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Abstract 

We study the two dimensional XY model with quenched random phases and its 
Coulomb gas formulation. A novel renormalization group (RG) method is devel- 
oped which allows to study perturbatively the glassy low temperature XY phase 
and the transition at which frozen topological defects (vortices) proliferate. This 
RG approach is constructed both from the replicated Coulomb gas and, equivalently 
without the use of replicas, using the probability distribution of the local disorder 
(random defect core energy). By taking into account the fusion of environments (i.e 
charge fusion in the replicated Coulomb gas) this distribution is shown to obey a 
Kolmogorov's type (KPP) non linear RG equation which admits travelling wave so- 
lutions and exhibits a freezing phenomenon analogous to glassy freezing in Derrida's 
random energy models. The resulting physical picture is that the distribution of lo- 
cal disorder becomes broad below a freezing temperature and that the transition 
is controlled by rare favorable regions for the defects, the density of which can be 
used as the new perturbative parameter. The determination of marginal directions 
at the disorder induced transition is shown to be related to the well studied front 
velocity selection problem in the KPP equation and the universality of the novel 
critical behaviour obtained here to the known universality of the corrections to the 
front velocity. Applications to other two dimensional problems are mentionned at 
the end. 



1 Introduction 



1.1 Overview 



Topological phase transitions in two dimensions, which are induced by the 
proliferation of topological defects, are naturally described via Coulomb gas 
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formulations. In the description of the Kostcrhtz-Thoulcss transition of the 
XY model, the corresponding Coulomb gas (CG) with integer charges is ob- 
tained as an effective theory for the topological defects of the model : the 
vortices [1,2]. Similarly, a collection of dislocations whose proliferation in- 
duce the melting transition in a two dimensional elastic lattice can be de- 
scribed by a Coulomb gas with vector charges which belong to the reciprocal 
lattice[3,4]. Finally, most of two dimensional statistical models, such as the 
Ising, Potts and Askin- Teller models, can be transformed into Coulomb gases 
[5] (see also the review [6]). In all these cases, the transitions can be studied 
by renormalization in considering the screened interaction between two test 
Coulomb charges. This interaction is logarithmic at large distance r, and the 
renormalization procedure consists in neglecting the higher terms in an ex- 
pansion in 1/r as being irrelevant [6,7]. This renormalization is thus valid for 
a dilute gas of charges, and is usually implemented by an expansion of the 
Coulomb gas partition function in powers of the charge fugacity y. This fu- 
gacity y — exp{—E,./T) is related to the core energy of the charges (i.e the 
local energy to create a defect or their chemical potential). Thus this Coulomb 
gas renormalization procedure is perturbatively controlled in the limit of large 
or equivalently at finite temperature in the limit of small fugacity. Sev- 
eral topological transitions in two dimensions have been successfully described 
using this CG technique [2,3,6]. 

Soon afterwards, several authors attempted to extend these techniques to mod- 
els with quenched disorder [8-10]. The randomness in the original statistical 
models translates into random fields in the Coulomb gas formulation. The av- 
eraged free energy (instead of the partition function) of this Coulomb gas is 
then expanded, as in the pure case, in powers of the fugacity y of the charges. 
Usual scaling techniques then allow to study the topological transitions in 
these disordered systems. 

It does not seem to have been realized at that time that these approaches rely 
on the crucial assumption of a uniform fugacity for the charges in the sample, 
or equivalently of a core energy spatially uniform over the sample. Although 
this assumption is natural for pure models, it is at least questionable in the 
presence of a random potential [11]. The randomness may favor the appear- 
ance of topological defects (the charges of the CG) in some sites where the core 
energy will be effectively lower than on other sites. At high temperature, the 
randomness of the core energy E^ is irrelevant as it is averaged out by ther- 
mal phase ffuctuations. In that case the large scale behaviour of the model 
can indeed be described by simply considering the averaged fugacity over the 
sample ^, as in [8-10]. However, as we show here, when the temperature is 
lowered, the spatial inhomogeneities in the local core energy become more 
and more relevant. As a result, approaches based on a single uniform fugacity 
are doomed to fail. A correct detailed description of the scaling behaviour of 
the site-dependent fugacities (or core energies) becomes then necessary to de- 
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termine the phase diagram and describe quantitatively the transitions in the 
random model. It is the purpose of this work to define a novel renormaliza- 
tion method which allows this description in a consistent and perturbatively 
controlled manner. We present here a detailed analysis, a shorter account of 
the method and results has appeared in [12]. Although we focus here on the 
random phase XY model, these techniques can be applied to a much wider 
class of systems which can be formulated as disordered Coulomb gases (e.g. 
see Section 7 for a random Sine-Gordon formulation). In particular, extensions 
to the melting of two dimensional crystals in presence of disorder are stud- 
ied in [13-15] and apphcations to the statistics of locahzed wave functions of 
electrons in random magnetic fields as well as entropic transitions in Liouville 
theory and other models are studied in [16]. 

1.2 Random phase XY model 

In this paper, we focus on the 2d XY model with randomness in the phase (see 
eq. (1)), originally studied in [9] in the context of XY magnets with random 
Dzyaloshinskii-Moriya interactions. The topological defects (vortices) of this 
model are represented by integer charges rij. at sites r. Two charges rij., rij.' 
interact for a large separation via the usual Coulomb interaction of strength 
J, with the corresponding energy — 2J 12^11^' ln(|r — r'\/a) where a is the typical 
size of the core region of the vortices, iethe short distance cut-off (as discussed 
in section 2). These charges also couple to a random potential K, which arises 
from the randomness of the XY model, via a term rirK which depends on the 
sign of the charge. The crucial property of this random potential is that it has 
long range logarithmic correlations: 

(K - K')^ = 4(7 J2 In |r - r'l + 0(1) 

the on-site variance being — A 2a In L where L is the system size. 

Rubinstein, Schraiman and Nelson [9] identified the most relevant vortices as 
the n — ±1 charges. They derived scahng equations for the stiffness J(Z), 
the disorder strength a{l) and a single fugacity for these charges y{l)- As in 
the pure case, the nature of the phase is determined by the behaviour of the 
fugacity : in a quasi-ordered phase, y{l) decreases with the scale while in a 
disordered phase it increases. In this last case, vortices appear at the scale ae'* 
at which their renormalized core energy Ec{l) is about zero, which corresponds 
to a fugacity of order 1 : y{l*) — 1 (see e.g [17]). Interestingly the divergence 
of the correlation length at criticality found in [9], ~ exp{cste/\T — 2), is 
identical to the result of Kosterlitz-Thouless for the pure system. The phase 
diagram was found to be reentrant at low temperature, and is shown as a 
dashed line in fig. 1 as a function of the renormalized value of a and J. 
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Fig. 1. Phase diagram of the two dimensional XY model with random phase fluctu- 
ations. J corresponds to the stiffness, T to the temperature and a to the strentgh 
of the randomness. The dashed line represents the transition line of Rubinstein, 
Schraiman and Nelson [9] between the quasi-ordered (XY) phase and a disordered 
phase. 

Recently, several authors [18,19] have proposed a modified phase diagram for 
the same model, where the reentrance of the phase transition of fig. 1 dis- 
appears. The main point, further developed in [20], is an energy argument 
for a single vortex at zero temperature in a finite size sample, in the spirit 
of the Kosterhtz-Thouless argument for the transition in the pure model [2]. 
The energy to create a single defect has two main components (neglecting 
the bare core energy) : the elastic energy Ef.i ~ JlnL and the disorder en- 
ergy. Indeed the defect can take advantage of the spatial variations of Vy. and 
choose the minimal value of K- in the N = (L/a)"^ sites of a sample of size 
L. To estimate this minimum Vmin oi Vr over N sites, the above authors 
neglected the spatial correlations of Vj.. Under that hypothesis, the problem 
of this single particle in the random potential Vr becomes identical to the 
Random Energy Model (REM) defined and solved in the seminal work of 
Derrida [21]. In particular, the averaged Vmin of the minimum behaves for 
large N as V^j„ ~ — \/2A(lnA^)2 ~ —^/SaJlnL, where A is the onsite vari- 
ance (see above) . Adding the elastic contribution E^i yields a creation energy 
E ~ J(l — \/8o") InL. Hence this simple single vortex argument points towards 
the existence of a topological phase transition atcr = o"c = = where 
vortices proliferate. The modified phase diagram is shown as a solid line on 
figure 1. 

Although appealing, this single vortex argument is not sufficient by itself to 
prove the existence of a phase transition in the system of many interacting 
charges, and even less so to describe the critical behaviour. Screening by mu- 
tual interactions usually plays an important role in Coulomb systems, and 
screening of the disorder could also in principle modify the results. The sin- 
gle vortex argument can only become valid asymptotically (at infinite scaling 
length) if there indeed exists a phase with finite renormalized values for J and 
a and vanishing fugacity at zero temperature. The very existence of the XY 
phase in the present problem has been questionned in recent works [22] and a 
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further analysis is thus necessary. Thus, we cannot avoid the use of a renormal- 
ization approach to determine the true phase diagram. Such a controlled RG 
description should correctly take into account the correlations of the random 
potential, neglected in the above argument. It should allow to precisely car- 
acterize the universal features of the critical behaviour at the disorder driven 
transition. 

Several attempts to construct a RG method were proposed [19,23-25] prior 
to this work, but no agreement was found between the results of these dif- 
ferent approaches. As we discuss below, this is largely because none of these 
approaches was internally consistent. In particular, they were all based on ex- 
pansions in dipole fugacities while considering single charge fugacities not only 
appears more natural but is essential in constructing a consistent renormaliza- 
tion procedure. As a result these approaches missed the very important contri- 
bution of the fusion of environments (see below). Korshunov and Nattermann 
started by approximating the free energy of a given disordered Coulomb gas by 
the sum of the free energies of independent dipole^ of charges [23] . They found 
no renormalization of the disorder strength a and a modified phase diagram in 
agreement with the figure 1. An interesting replica approach was developed by 
Scheidl [24], with the use of a Coulomb gas of m-component charges n, where 
m is the number of replica. Scheidl noticed that, in the renormalization proce- 
dure, one must a priori take into account charges with a number p > 1 of non 
zero components, contrarily to [9] where only single component charges were 
considered. His RG equations are not compatible with the one of [23]: since 
the disorder strength a is renormalized, it leads to a different phase diagram 
when expressed in the bare constants. Upon closer examination in Section 8 
and interpretation within our formalism, the assumption implicitly underlying 
Scheidl's work [24] appears to be the gaussian nature of the renormalized local 
disorder, while we are able to prove that the local disorder does not remain 
gaussian upon coarse-graining. Although we are also able to show a posteriori 
that in the XY phase some results are compatible between the two approaches, 
the non gaussian nature of the local disorder becomes crucial to determine the 
critical behaviour at the transition. Finally Tang [25] developed an appealing 
physical picture for the single vortex freezing phenomena in this model be- 
yond the REM approximation mentionned above, and correctly foresaw the 
existence of a connection [25,26] with the problem of the Directed Polymer on 
the Cayley Tree (DPCT). However, no precise Coulomb gas renormalization 
procedure was developed following these ideas. 

Finding a controlled RG procedure which allows to describe this physics of 
freezing in a collection of interacting Coulomb charges in a random environ- 
ments is thus a non trivial and challenging problem. As we will see (and as 
can be expected from the above single charge argument), it requires a quan- 
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titative description of the scale dependence of the probabihty distribution of 
the fugacities associated with the rare sites were the charges are frozen. 

Before embarking in the (sometimes technical) remainder of the paper, it is 
useful to depict the spirit of the method and summarize the main results. 



1.3 Description of the methods 

To capture the physics of freezing in a gas of Coulomb charges we find it 
crucial to start by realizing that upon increasing the size of the charges (cut- 
off a), their core energy aquires a random component from the potential Vr. 
Increasing a to a = a(l + d/), this potential K with logarithmic correlation 
splits into the rescaled logarithmically correlated potential at scale a, and 
a random part Vj. uncorrelated at scales larger than a. This local potential 
is naturally incorporated in the (rcnormalized) core energy Ec, which thus 
becomes random and site dependent : E^'' E^"^ = Ejf"^ ± Vy. Thus, upon 
coarse graining, the fugacities of the ±1 charges of the Coulomb gas become 
random (and dependent on the sign of the charge) : z± = ye^^'"". Let us stress 
that the definition of Vr (and thus of z±) depends on details of the cutoff 
procedure. However the RG procedure developed here depends only on the 
logarithmic behaviour of the correlator of at large distances, which is cut- 
off independent. As a result, a remarkable universality will emerge at the end 
of the procedure, which happens to be related to universality of front solutions 
in non- linear equations (see below). 

We are now faced with the description of the scale-dependence of a Coulomb 
gas with random fugacities, and thus a priori of the full corresponding fugacity 
distribution. Since we find that, at low temperature, the distribution of the 
local core energy does not remain gaussian under coarse graining we cannot 
restrict a priori the renormalization study to follow a small number of variables 
(such as the mean and the variance) as illustrated in figure 2. Instead, we must 
study the scale dependence of the full distribution P{Ec) and especially the 
precise form of its tails around < (i.e |/ ~ 1) which control the low 
temperature physics (at T = defects appear where the local core energy is 
negative). This requires new techniques in two dimensions. 

Before developing a new RG procedure, we have to find a correct perturba- 
tive parameter to study the disorder driven transitions at low temperature. 
Around the pure topological transition, this parameter corresponds to the 
charge fugacity y (for the most relevant charges n = ±1). However, we expect 
that, upon lowering the temperature, a freezing of the defects occurs. In that 
case, most sites have a large core energy {y ~ 0) while only a few sites with 
Ec are favorable to the defects. A natural choice for a perturbative pa- 
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Fig. 2. Schematic representation of the probabihty distribution P{y) of the local 
fugacities y, plotted as a function of — In y = /5£'c (where Ec is the renormalized 
core energy). At high temperature (left) it is narrow and a description in terms of 
the averaged fugacity is correct. At low temperature the distribution of /SEc is broad 
and non gaussian and information about the whole distribution is a priori needed 

rameter at low temperature in this highly non-homogeneous Coulomb gas is 
thus the density of favorable szte^ which we note P{y ~ 1). Pi{y ~ 1) plays 
an analogous role here than y{l) in the pure case: the quasi-ordered phase 
is caracterized by a density Pi{y ~ 1) decreasing with the scale /, while a 
disordered phase, where disorder induced topological defects prohferate, cor- 
responds to an increasing Pi{y ~ 1). This new small parameter allows us to 
study in a controlled perturbative expansion the physics of the vortices at low 
temperature and around the transitions. 

An essential contribution to the renormalization of Pi{z^,Z-), absent in pre- 
vious approaches, originates from what we call the fusion of random environ- 
ments. The main idea is the followoing. Each random fugacities is associated 
with a region of size a around r. Upon increasing the cut-off, the size of these 
regions increases, and we thus have to merge two regions distant from less than 
the new cutoff a (see fig. 3). The probability distribution of the fugacities in 
the new region can be deduced from the one in the two merged regions by a 
fusion rule (see fig. 3). The final renormalization of Pi{z+, zJ) can be formu- 
lated into a single differential equation for the distribution function P;(-2+, z_) 
(given in Section 3). Most interestingly we find that the universal features at 
the transition of the model (a^ = |) do not depend on the precise definition 
of these regions. The fusion of environments corresponds to a non-linear term 
in the differential equation, which can thus be affected by a change of cut-off 
procedure. Quite remarkably, the universality class at the new transition is de- 
termined by the properties of the front solutions of this differential equation 
(i.e their velocity), that do not depend on the precise form of this non-linear 
term. Hence universality in the usual RG sense finally appears from non trivial 

^ Actually, the random fugacities y^ depends on the sign of the charge : z^-, and 
the perturbative parameter is the density P(z+ ~ 1) = P(z_ ~ 1) of sites favorable 
either to +1 charges or to —1 charges. One must thus follow the RG flow of the full 
probability distribution P;(2;_|_, z_). 
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Fig. 3. Schematic representation of the fusion of two local environments upon 
coarse-graining. Two regions, distant from less than the new cut-off d, arc merged 
when the cut-olf is increased. The corresponding fusion rule for the distribution 
function of the local variables is shown on the figure. 

results of front propagation in non-linear equations. 

This differential RG equation for P/(-2+, z^), together with the screening equa- 
tions for J and a, can be obtained by a systematic and perturbatively con- 
trolled RG procedure. One formulation consists in using the replica trick to 
average the free energy over disorder, leading to a Coulomb gas with m- 
component charges. In this formulation, the infinite set of fugacities, associated 
with the vector charges Y[n], encode the distribution function Pi(z+,Z-) of 
the random fugacities z^. Hence, since this last distribution cannot be a pri- 
ori described by a finite number of moments at low temperature, we have 
to consider the scaling behaviour of the fugacities of all the replica charges 
with components 0, ±1. The m — > limit of this infinite set of RG equations 
is then taken by using an appropriate parametrization in terms of Pi{z^,Z-) 
which yields the non linear RG equation for -2-). The second, equivalent 

formulation, does not rely on replica. It is constructed using a new expansion 
of physical quantities in the "number of independent regions" , introduced in 
this paper. 

1.4 Disorder induced topological transitions and connection with DP 

After some transformations, we reduce the non-linear RG equation that gov- 
erns the scale dependence of the distribution Pi{z^,z_) to the celebrated 
Kolmogorov-Petrovskii-Piscounov (KPP) equation (also named the Fisher equa- 
tion). The known results on this equation allow us to derive the form of the 
tails of the distribution and thus to obtain the phase diagram. In particular, 
the velocity of the front solutions of the KPP equation directly determines 
whether the small parameter Pi{y ~ 1) increases or decreases. Hence using 
known results on the selection of this velocity, we obtain the phase diagram of 
the figure 1 expressed in renormalized variables. The critical behaviour at the 
transition cr = |, follows from the finite size corrections to the velocity cor- 
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rections of the KPP front and defines a new universality class. In particular, 
the correlation length ^ diverges at the a — ac — l,T — transition as 



in contrast with the KT behaviour found in [9]. Note that recent numerical 
simulations [27,28] of this model seem to agree with the phase diagram of fig. 
1 and it would be interesting to also determine numerically the precise critical 
behaviour at low temperature. 

Besides the caracterization of the new critical behaviour our work also en- 
ables to study the freezing of vortices, which occurs below a temperature Tg 
(see figure 1). This freezing corresponds to a transition for the single charge 
problem, whose study with our new RG technique is presented in [16]. In par- 
ticular, this renormalization approach draws a precise connection between the 
physics of freezing of the XY defects and the problem of branched processes 
(or of random directed polymers on Cayley trees (DPCT) for a discrete ver- 
sion) studied by Derrida and Spohn [29]. This connection naturally emerges 
from our Coulomb gas RG equations via the KPP equation which has also 
appeared in the exact solution of the DPCT problem [29] . It does not rely on 
any ad-hoc construction. 

The paper is organized as follows : in Section 2, the random XY model is 
defined and its CG formulation is carefully derived. In particular, the relation 
between the continuum limit and the decomposition of the random poten- 
tial 14 in to a local part (random core energy) and a long-range potential is 
discussed in section 2.2. Part 3 describes the renormalization method of the 
replicated Coulomb gas, and while a direct method (without replica), which 
consists in expanding the free energy into the number of independent regions 
(random environments) is presented in Section 4. The RG equations are an- 
alyzed in Section 5 using results on the propagation of KPP-like fronts. The 
consequences for the determination of the phase diagram and the critical be- 
haviour is detailed in Section 6. A formulation in terms of a Sine-Gordon model 
is given in Section 7. The comparison with previous approaches is postponed 
to Section 8, and most of the technical details can be found in the appendices. 



2 XY Model with random phases 
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2.1 lattice Coulomb gas 



In this Section and in the following we study the XY model with random 
phases [9]. We start with the model defined on the 2D square lattice by its 
partition function: 



z[A] = n / de, e-^^i'^'^i with H[e, A] = J2 v{e. - e, - a,) (i) 

where the sum is over pairs of nearest neighbors (i.e over bonds) on the lattice 
and (3 = 1/T is the inverse temperature. The Aij are random gauge fields, 
independent from bond to bond, each with gaussian distribution of variance 
Afj = 71(7. The periodic potential V{9) is defined for the XY model by V{9) = 
— - cos(6') where J is the stiffness. In the limit a ^ oo this model corresponds 
to the 2D "gauge glass model" [30,31]. For finite a it was studied in Ref. 
[9,25,24,23]. 

The standard way to study this model is to decompose it into spin waves and 
vortex degrees of freedom Z[A] = ZswZcg- This decomposition can be per- 
formed exactly (see Appendix A) for the corresponding Villain model defined 
by the potential 

+°° 

p=— oo 

Technically this is the model which we study here. It is reasonable however 
to expect that this Villain model and the XY model should be in the same 
universality class (as is the case without disorder). The RG analysis contained 
in the following Sections is consistent with this assumption . The vortex 
part is described in terms of a Coulomb gas with integer charges rij. defined 
on the sites r of the dual infinite (square) lattice: 



ZcG = ^latt = E e"^"" (2a) 

{Ur} 

H = E ^Jrij.Gr-r'IT'r' — E ^rK (2b) 

2 r^r' r 

where Gr-r' = /kG'k(l — e*'' *^'"^'"'^) is the 2D Coulomb potential with Gr-r' ~ 
In |r — r'l at large distance. G^^ = ^[2 — cos{kxao) — cos{kyao}] is the lattice 

Laplacian and we denote = Jli"/°a^ J'-i%o ^{2nf where Oo is the lattice spac- 
ing. Note that the energy associated with a dipole of unit charge of size r is 

^ It can be shown to lead to unimportant additional random terms in the bare 
fugacity 
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2Jlnr. As discussed in the Appendix A we can consider a neutral CG (since 
J q = 0) with J2r^r = 0, and neutrahty has aheady been used to arrive at (2). 

In the vortex representation the random gauge fields A translate into random 
dipoles (along z) qr' = ^^r' x A which couple to the vortex charges via 
the Coulomb potential. As detailed in Appendix A this results in a gaussian 
bare disorder potential Vj. = — 2JX]r' ^r-r'Q'r'- An important feature of this 
problem is that the disorder potential seen by the charges (the vortices) has 
logarithmic long range correlations 



(Vr-K')^ =4ajMn|r-r'| +0(1) (3) 

for r 7^ r', since VkVCk = 2aJ'^G\^. Note that at a given point ^ 2aJ'^\n.L 
where L is the system size and that up to now, these are exact transformations. 

Before defining the continuum limit of this Coulomb gas, we note that an 
alternative definition to this lattice Coulomb gas consists in labelling the non 
zero charges in (2) by their positions rj and their corresponding charge n,. 
Instead of the integer field Uy. of (2) defined at each site of the lattice, a con- 
figuration is represented by the set {n,, rj}. We obviously have Uy. = J2i ^A,ri- 
With this representation, the partition function of (2) reads 

p>0 {ni,...np} ri^---7^rp 

where Gr^-rj and Vr- have been defined after eq. (2) and the sum over the 
charge configurations (primed sum) counts each distinct neutral configuration 
of non zero charges only once[^. 



2.2 continuum limit and decomposition of the disorder 



In order to implement a renormalization procedure one first needs to introduce 
a continuum version of this model. In the usual approach to 2D Coulomb gas 
the continuum limit is obtained by replacing the lattice Coulomb interaction 
by the approximation [2,32] 




Gr-r' ^ G:Z' = In ^^-^ +7(1- ^;^) (5) 



This leads to the factor 1/ N{n)l in the definition of the configuration sum of 
the CG, where N{n) is the total number of particles of charge n 
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where 5^"°/ = 1 for |r — r'| < Uo and otherwise, and 7 = ln(2-\/2e'^) with 
C = 0.577216 is the Euler constant. This approximation is excellent [33] on 
the lattice for |r — r'| > Oq. In the standard method the continuum CG model 
is then defined by considering a gas of integer hard core charges Ur- at point rj 
of diameter Ug which interact with G^^^,. Using neutrality X^r^r — and (5), 
allows to rewrite [2,32] the hamiltonian (4) as a sum of a simple logarithmic 
interaction —J Y^ij^j ^r, In ^'"'"'"^ ^ ) rirj between the hard core charges and a 
fugacity term InyX^jn^. of bare value y = e'^^" where = 7 J can be 
interpreted as the bare core energy for the defects of the model (2). 

In presence of disorder, special care has to be taken to define properly the 
continuum limit for the random potential, since its correlator is logarithmic. 
Since on the lattice the correlator of the disorder G is the same as the Coulomb 
interaction, it is consistent to use the same G^pp for the disorder in the con- 
tinuum model. This immediately leads to the fact that the disorder K must 
be separated in two parts, using (5): a long range correlated gaussian part 
and a local part Vj-: 



V, = V> + V, (6a) 
{V> - V>Y = AaJ' In f (1 - 5^) (6b) 



= 2(7 727^";^ (6c) 

with no cross correlation. Using this decomposition we can now write the 
partition function of the continuum model: 



-..-l+E E / !^...!^e-*''M (7) 



'-_P^-pH[n,r] 



H[n, r] = - J ^ n. In ^ ^ - ^ n,V> - ^ In F[n„ rJ (8) 

i^j \ J i i 

where the primed configuration sum counts only once each distinct neutral 
charge configuration. We have introduced the spatially dependent fugacity, of 
bare value, from (5,6): 



In y[n, r] = --fpJn'^ + pnv^ (9) 

Thus we find that disorder favors some regions resulting in a local fugacity 
for ±1 charges t/(r) = e^^^"^^^ with a core energy -E'c(r) = E^^Vr which now 
varies from point to point. Thus one anticipates that problems will arise in 
the conventional fugacity expansion if |/(r) varies strongly from point to point. 
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In addition, note that the local fugacitics are different for + charges and — 
charges (although there is still a statistical +|— symmetry). 

We have thus defined a continuum model with a particular cutoff procedure. 
This is apparent, e.g. in the decomposition of the disorder that we have used. 
This decomposition is cutoff dependent and we have taken care to chose the 
same cutoff procedure for the disorder term and the interaction. We have 
chosen here a real space hard cutoff procedure, which is often used in CG 
studies. Other cutoff procedures can be used. It is natural to expect, and 
we will partially verify that the large scale results will not depend on the 
particular procedure chosen. 

The alert reader will have already noticed that the disorders and Vr as 
defined above in (6) are not strictly speaking physical, since the Fourier trans- 
form of their respective correlators is not positive. This is an artefact of this 
particular choice of a real space hard cutoff. It is not a serious problem, and 
is easily cured by choosing instead a cutoff in Fourier space. The resulting de- 
composition of disorder is then completely legitimate. This is further explained 
in Appendix B. For simplicity, we will however proceed using the above cutoff 
choice, which has illustrative value, keeping in mind that the more legitimate 
choice detailed in Appendix B can be used instead, completely equivalently at 
all stages, for technical rigor. 

We now turn to the renormalization of the model. 



3 Renormalization using replica 

In this Section wc study, using replica, the renormalization group properties 
of the disordered Coulomb gas defined by (2). In the present case the replica 
method is particularly convenient in order to perform the combinatorics nec- 
essary to renormalize consistently the model. 

The strategy is first in (3.1) to transform the model (2) into a vector Coulomb 
gas with m-replica charges. The fugacities of these m- vector charges will then 
naturally encode the distribution of the spatially dependent fugacities defined 
above (see (9)). The renormalization group (RG) equations for these fugacities 
are derived in (3.2) for fixed m. By a suitable parametrization, in (3.3) we 
then extract in the m — > limit the RG equations which describe the scale 
dependence of the full distribution P{z+, zJ) of the local fugacities of the 
topological defects in the original disordered model (2), as well as the scale 
dependence of the stiffness J and long wavelength disorder a. 
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3. 1 Replica Coulomb gas, continuum limit and fugacity expansion 



We start again from the model (2) on a square lattice. As is well known, 
disordered averaged correlation functions and free energy can be obtained by 
studying the replicated partition function (generating function) Zl^^t in the 
limit m — s> 0. For integer m, Z^^a can then be written exactly as a CG with 
m-component vector charges n^, a = 1, ..m (each is integer), living on the 
sites of the lattice (note that a denotes a replica index while "co" denotes 
the cutoff). Averaging over the bare disorder one obtains the partition sum 
of a fully coupled, translationally invariant, vector Coulomb gas on a square 
lattice: Z^^, = E{ng} e'^^i-t with 

= - E <a-r'n^ (10) 

r^r' 

where Kab = l3J5ab — o-p^J^ and summation over repeated replica indices is 
assumed unless otherwise specified. 

The next step is to approximate the lattice replica model (10) by a continuum 
Coulomb gas with m-component vector charges. In the following we consider 
a hard core cutoff in real space. The problem of the choice and consequences 
of the cutoff procedure is rather subtle here and will be discussed below. Using 
the approximate propagator (5) we obtain the continuum hamiltonian 



r] = - ^ KabntGt^l^M'j (11a) 
= - 5: In (h^] n] - ^ In F [n,] (lib) 



where the charge rij is located in r^. In the second equality wc have used the 
neutrality of the Coulomb gas, i.e. X^r^r — f*^^ each a = 1, ..m, to introduce 
the local fugacity Y[\\\ which is a function of the whole set of components of 
the vector charge n = [n^, ..rf^). Its bare value from (5) is a simple quadratic 
function [24]: 



l^[n]bare = e"""^^"'"^ (12) 

This quadratic form results from the Gaussian nature of the hare local disorder 
and corresponds to (9) in the unrcplicated version. If this form was preserved 
by the RG, as was implicitly assumed in [24], one would be able to study the 
model using only two coupling constants. However this is not the case. As 
shown below, the vector charge fugacity y [n] has a non trivial flow under RG 
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and does not remain purely quadratic. The local disorder does not remain 
gaussian and we will have to follow its full probability distribution. 

We now study the scale dependent properties of the m-component vector 
Coulomb gas using the expansion in the vector charge fugacity Y[n]. Although 
it is the natural way to study the renormalization of a vector Coulomb gas, it 
may seem at this stage somewhat formal. This is not so however since, as will 
become clear below, it turns out to correspond exactly to the expansion in 
the number of rare favorable regions of local disorder, which is the physically 
relevant (and novel) expansion for this model. For the replicated partition 
function in the continuum model this vector fugacity expansion reads: 



EE/,, '-^...'-^e-^"'^^-" (13) 



which contains fugacities via (11) and the primed sum over charge configu- 
rations counts each distinct neutral configuration only once. In presence of 
disorder, infrared divergences appear everywhere in the low temperature XY 
phase [34] . To treat these divergences we now turn to the RG method. 



3.2 RG equations 



We perform the RG analysis of the present m-component vector Coulomb gas 
on the partition function Z"*. It is a simple extension of the analysis for the 
scalar Coulomb gas [2,6]. Details arc presented in the Appendix C and we 
only sketch the method here. For any fixed m it is possible to leave the form 
of the expansion (13) unchanged under the increase of the hard core cut-off 
Oo OoC'^^ provided one defines scale dependent coupling constants K^'^ and 
fugacities Yi[n]. This corresponds to the (one loop) renormalizabihty of the 
m-component vector model, which we checked here to order y[n]^. The RG 
flow equations which determine these couplings are found as (see Appendix 
C): 



di{Kf%, ^c^J2 nVr[n]r[-n] (14a) 
117^0 

diY[n] = {2-n''K,hn')Y[n]+C2 ^ Y[n']Y[n"] (14b) 

n'+n" = n 

with Ci = 27r^, C2 = vr for our hard cut-off procedure, and the second equation 
(14b) is defined only for n 7^ 0. The first equation (14a) comes from the 
annihilation of dipoles of opposite replica vector charges separated by Go < 
gives the renormalization of the interaction (screening by 
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Fig. 4. Fusion of replicated charges. In cases where the only charges which are 
relevant are the single component ±1 charges the fusion (a) always leads to charges 
with a higher number of components ±1 and are thus irrelevant: in that case the 
fusion upon coarse graining does not play any role. On the contrary, in the present 
case multicomponent charges are also relevant and must be taken into account: 
fusion of charges (b) does not necessarily increase the number of ±1 components 
and must be taken into account in the renormalization of the fugacities of the 
relevant charges, (c) (and (a)) illustrate examples of fusion of vector charges in the 
replicated Coulomb gas which do not correspond a true fusion of scalar charges in a 
single environment (nor to any annihilation as in (b)). Instead they correspond to 
the "fusion of environments" and encode for the renormalization of the probability 
distribution of local fugacities upon coarse graining (as they keep track of some non 
trivial correlations, see discussion in the text). 

small dipoles) and of the disorder. Simple rescaling gives the first term of the 
second equation (14b), i.e. the naive scaling dimension of Y[n]. 

The second contribution in (14b) comes from the possibility of fusion of two 
replica vector charges upon coarse graining (see fig. 4). While such term can be 
neglected in the (pure) scalar Coulomb gas (as it yields less relevant operators) 
it is usually crucial when studying most vector CG models, as e.g. in the 
analysis of two dimensional melting transition [3] . Indeed following too closely 
the analysis for the pure XY model has led previous studies [9,19,25,24,23] 
to miss the possibility of fusion and thus such a contribution. As we see in 
the following it has important and non trivial consequences for the physics 
of the low T phase and the transition. For the present disordered Coulomb 
gas, contrarily to the conventional analysis [3], one cannot hope to capture 
the most relevant operators by restricting to single component charges (e.g. 
n"" = ±()aai)- This was recently emphasized by Schcidl [24]. However, since 
this leads to considering multicomponent vector charges, it is thus crucial to 
treat properly this fusion term, which was not done previously (e.g. in [24]). 
Moreover, discarding this term in (14b) leads to a set of RG equations which is 
not consistent to their lowest order 0(y[n]^). This term may a priori modify 
the scaling dimensions in a non trivial way in the m — > limit, and it is thus 
crucial to study carefully its effect. 

Before doing so in the next Section, let us give for completeness the renormal- 
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ization of the free energy density per replica, defined as fm — a^L InZ"^. 
It reads: 

from which the flow of the free energy density can be obtained as f — 

l™m-»0 fm- 

3.3 limit m — >• and fusion- diffusion formalism 

We now have to find an analytical continuation to m ^ of the whole set 
of RG equations (14). This is a priori a formidable task because (14) are in 
fact, for arbitrary m, an infinite set of coupled equations. Remarkably, in the 
process of performing this analytical continuation, an appealing physical inter- 
pretation in terms of probability distributions of local fugacities (local disorder 
: see Section 2.2) will emerge naturally and be our guide in the following. 

3.3.1 parametrization of the fugacities and analytical continuation 

As a first step, wc will consider only charges with components = 0, ±1 in 
each replica. How to incorporate higher charges (e.g. = ±2 . . . ) is discussed 
in Appendix 1, where it is shown that they are less relevant in the region of 
the phase diagram studied here. We first remark that the possible forms of the 
Y[n] are severely constrained. Replica permutation symmetry, which we will 
assume here and is preserved by the RG, together with ~ 0, ±1 implies that 
Y[n] depends only on the number n+ and n_ of +1/ — 1 components of n. A 
natural possible parametrisation of Y[n] = y[n+,n_] consists in introducing 
a function of two arguments $(-2+, Z-) such that: 

y[n] = {z+'^z1~)^ = (JJ [Sno.,0 + Z+Sna,+l + zSna-l])^ (16) 

a 

where we denote {A)^(^z) = Jz+ z- ^ ^i^+j ^-)- strategy is to establish an 
RG equation for ^{z+,z_) (in the limit m ^ 0) whose solutions ^i{z^,z_) 
will parametrize solutions Yi[n] of (14a, 14b). 

Let us now examine how (14b) can be transformed in an intcgro-differential 
equation for $(z_|_,z_). The technical details are given in Appendix D . The 
first terms in the r.h.s. of (14b) translate into a differential operator (2 + (9)$ 
where: 
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O = /3J(2 + z+d,^ + z^d,_) + aipjfiz+d,^ - z^d,_f 



(17) 



Using new "core energy" variables u^v such that z± = e^^'^'^'"\ and the corre- 
sponding function $(m, v) such that $(it, v)dudv — $(-2+, -2_)(i-2+(i2;_, it can be 
interpreted as a diffusion process since the first term of (14b) now translates 
into (2 + (5)l> with O = ja„ + aJ'^dl. 

To deal with the second term we first extend the RG equation (14b) so as 
to allow for zero charge n = 0, since it is easier to continue analytically 
unrestricted sums. After some combinatorics (sec appendix D), we find that 
using the representation (16), (14b) can be rewritten completely equivalently 
in terms of $ as: 



di^{z+, zJ) = (2 + 0)^{z+, z_) - 2c2U^z+, z_) + C2NH{z+)5{z_) (18) 



/ z' + z" \ ( z'_ + z"_ 



1 + z'_z\ + z'^z'L ) \ 1 + z'_z\ + z'^z'L ) 

where M = J^^ ^_ $(-2+, z_). This equation describes the scale dependence (/) 
of the function ^i{z+,z_) which parametrizes the whole set of scale dependent 
fugacities Yi[n] in the limit m — > 0. 



3.3.2 RG equations and fugacity distribution 

Up to now the function $(2;+, z^) has been introduced as a generating func- 
tion to parametrize the fugacities Y[n]. It is a priori an arbitrary function 
and in particular J\f = / $ is still undetermined. In this paragraph we will 
exchange ^{z+,z_) for a physical function P{z+,z_) of norm unity, which 
will be interpreted in the following as the probability distribution for the lo- 
cal fugacities z^,z^ of ±1 charges. We start from the above equation (18) 
for ^(z-^-,Z-) which can be simply interpreted as describing the sum of two 
processes. Defining from the random fugacities z± = e^'^"^''^) the random core 
energy variables Ef = —{u±v), the first process in (18) corresponds to a 
brownian diffusion for v (i.e the local disorder potential as in (9)) together 
with a convection for u. The second process involves a fusion, with a rate C2 
upon increase of the cutoff, of two sets of random variables (z'^, z'_), {z'^, z'!_) 
into a single one (2;+, Z-) according to the transformation law: 

as in a A + A — A reaction. The term — 2c2A/'$ in (18) corresponds to a 
loss of two charges, while the last term corresponds to a creation of the fused 



18 



one. The term (){z^)S{z+) keeps track of the "dead charges" which disappear 
by setting them to (since they decouple from the system). It is in a sense 
only a counting device, since by construction J\f = / $ is unchanged upon 
fusion. We thus introduce $>(2;+,;2_) restricted to z+ > 0, > 0, such that 
$ = + (A^ - Af^)S{z-)S{z+) where Af> = J,^ 

>o,z_>o^ is total weight 
of non zero charges. Integrating (18) over z+ > 0, Z- > we obtain that: 

diX> = 2M> - C2MI (20) 

Thus in the presence of fusion it converges quickly towards — 2/ 02- 

Since H> = / $> converges to a constant, this suggests to introduce a normal- 
ized function P{z+,Z-) — $>/j\/'>. As shown below it is natural to interpret 
P{z-\-^ Z—^ as a probability distribution. Prom (18) we find that it obeys the 
following RG equation: 



diP{z+,z^)^OP-2P{z+,z^) 

where {..)p'p" denotes J^,^^^, ^//^^// P{z'_^_, z'_)P{z'[, z'i) and the probabihty con- 
serving diffusion operator O has been defined in (17). 

The limit m — * of the other RG equations (14a) which give the renormal- 
ization of the stiffness J and disorder strength a is performed in Appendix D 
using Reexpressed in terms of P they read: 



dj-^ _ 4ci / z'^z'L + z'_zl + ^z'^z'Lz'_zl \ .^^ , 



dl 4 \ {1 + z'+z'l + z'_z'^) 



pp 



da _ Ac, I {z'^z"_-z'_z\f \ ^22b) 



dl cl \{l + z'+z'L + z'_z'if 



pp 



where we have chosen (arbitrarily) to keep T fixed and renormalize J (only 
the combination K — (3 J flows). 



The above formulae (21,22) forms our complete set of RG equations. As will 
become clear in the following Sections, P(2;+,2;_) represents the distribution 
of the fugacities 2;+, z_ of local environments and the last term in (21) cor- 
responds to fusion of environments upon coarse graining. Remarkably, once 
expressed in terms of P the coefficients of the above RG equations exhibit some 
universality. The factor of 2 in the last term in (21) arises from the fraction of 
environments diV/V = 2 which are fused when increasing the cutoff. Note also 
that the coefficient ci/c^ which naturally appears in (22) is not affected by a 
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uniform rescaling of the fugacities. Note that some features of these equations 
are cutoff dependent, as will be discussed in a following Section. 

Finally, the flow of the free energy density is found to be: 



These RG equations will be studied in Section 5. First we will present another 
renormalisation procedure, without replicas. Although it is technically more 
difficult to implement, it allows for a more direct physical interpretation, which 
in turns sheds some light on the more systematic replica method presented 
above. In addition it may be more appealing to the replicaphobic Section of 
the community. 

3.4 validity of the method and universality 

The above method which relies on an expansion in the vector fugacities Y[n] 
can be justified provided there are few vector charges in the system i.e in the 
dilute limit. Even though the Y[n\ may appear as formal fugacities, the above 
RG equations can be justified in an expansion of the exact renormalised poten- 
tial G^{r) seen by two test charges distant of r in ao/r, as was emphasized by 
Nienhuis [6]. It is even claimed to be exact [6] in that limit[^. As will become 
apparent in the following Section the physical meaning of this diluted limit of 
vector charges exactly corresponds to the limit of a small density of regions 
favorable to the creation of frozen defects which is the physically relevant limit 
in the regimes studied in this paper. 

To understand how cutoff dependence comes in the method used here, it is 
instructive to study the limit of zero disorder. One can indeed check that one 
recover the usual results in the limit of the pure case a = 0. This however, 
requires some careful consideration of the cutoff procedure for the vector CG 
representation. As discussed in the Appendix F, even in the pure case the 
distribution $(-2+, Z-.) which parametrizes the vector fugacities solution of the 
CG RG equations can be non trivial. It does satisfy P{z+, z^) ~ S{z+ — z_) 
but the fugacity Zj^ = Z- = y still has a non trivial "distribution" $(?/) for a 
generic choice of cutoff. There is no paradox there and it is compatible with the 

^ In that case a cubic term should be added to them, — /iy[n] X^n/^o ^['^']^[~'^'] 
with h = vr(|7r + \/3)/2 in the second equation. We omitted this term as we found 
that it vanishes in the limit m ^ 



dif = 2/ - Tvr^ (ln(l + z'^z'L + z'_z'[)) 




(23) 
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standard Korsterlitz-Thouless RG equation of the pure case, as is explained 
in Appendix F, universality being recovered at small fugacity. 



4 Direct method of renormalisation 

In this Section we introduce a method to study the model (2) and more gen- 
erally Coulomb gas with disorder without using replicas. We start in Section 
4.1 from the general motivation and rederive the RG equation for P{z+,Z-) 
in (21) in a more physical way. We also identify the small parameter which 
allows to study perturbatively the present problem. In Section 4.2 we intro- 
duce a quantitative and systematic method to expand in this small parameter. 
The direct renormalization approach using this expansion is performed in Sec- 
tion 4.3. The connection with the replica method of Section 3 is presented in 
Section 4.4. 



4-1 RG method for broad disorder: physical derivation 

Let us first explain the spirit of the direct method and illustrate how one is led 

to the RG equation (21), derived more quantitatively in the next Sections. We 
have seen in Section 2.2 that the local disorder Vr defines the site dependent 
fugacities. We concentrate on ±1 charges for which these fugacity variables 
read (see (9)) 

zf = yrexp{±(3vr) (24) 

and are quenched random variables with only short range spatial correlations. 
One now studies the system under a change of cutoff ao aoe*^' (coarse grain- 
ing) which includes an integration over the corresponding degrees of freedom. 
We find that the coarse grained model remains of the same form as the orig- 
inal one, with a renormalized stiffness J;, a renormalized gaussian long range 
disorder strength ct; Jf and a local disorder distribution Pi{z+, z^). Note that, 
although the bare local disorder f r is gaussian, it becomes non gaussian under 
coarse graining. This is a novel feature of the present approach, at variance 
with previous attempts at renormalizing the model [23,25,24]. It complicates 
the analysis but is necessary to capture correctly the physics of the model 
which is driven by the rare events. 

The RG equations (21, 22) for the fugacity distribution Pi{z^,Z-) of local 
environments (higher charges are less important and considered later), for 
the stifness J; and for the correlated disorder strength ct; can be understood 
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from the following considerations. The correction to the fugacity distribution 
Pi{z+,z^) is the sum of two contributions: 

(i) rescaling: the first observation is that upon changing the cutoff, as can be 
seen from its correlator (6) and is detailed below, the long range disorder 
produces an additional local disorder contribution which can be written as a 
renormalisation of the local charge fugacity: 



^ z^e^i-^^^d-r^ (25) 

where dv,. is a gaussian random variable, uncorrelated from site to site and 
with dVrdVr' — 2J'^ This contribution leads to an effective diffusion 

and drift in the random core energy variables Ef — —{u±v) — —Tin 2;=*= as 
diEf — Jdl ^ dv and thus produces the first terms in (21). 

(ii) fusion of environments. The second contribution comes from the fusion 
of environments. Upon a change of cutoff, any two regions located around ri 
and r2 with oq < |ri — < aoc^' have to be considered as a single region in 
the system with the rescaled cutoff. As a consequence the two corresponding 
pairs of fugacities zf^ and zf^ must be combined and replaced by a single 
pair of effective fugacity variables associated with the new region at 

r = |(i"i + ^2)-, as illustrated in Fig. 5. z^ can be determined by estimating 
the relative Boltzman weight W+/Wo to have a configuration with charge 1 
(which lies either in ri or r2) versus a neutral one (either no charges or a 
dipole in ri,r2), and similarly for z_. This gives the fusion rule: 



= '"^ + (26) 

The corresponding correction to the distribution Pi{z+, z^) produces the last 
two terms of (21). 

Finally, the RG equation for J; and ai can be obtained from the screening by 
small dipoles of the effective interaction and disorder between two infinitesimal 
test charges as described in Section 4.3.3. 

Several comments are in order concerning this RG procedure. First we note 
that in defining local fugacity variables (24) we have added an explicit spatial 
dependence to the part i/r of the fugacity which does not distinguish between 
a +1 and a —1 charge. This dependence is not explicitly present in the bare 
model formula (9) (although it is present if an additional small disorder in 
the local stiffness Jij is included) but, as we can see from (26) it appears as 
soon as fusion takes place (the fusion rule is not compatible with a uniform 
Ur — y). Second, there are some assumptions underlying the RG procedure: 
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Fig. 5. schematic representation of the fusion of two local environments: when the 
cutoff is increased the two regions centered around ri (with local fugacities ^;!_) 
and r2(with local fugacities z", z") with qq < |ri — < aoe' are fused into a single 
region around r = ^(ri + with fugacity z+,Z- 

technically we treat the local regions as independent from point to point, we 
restrict V^-* to be strictly gaussian, together with the usual assumptions (e.g. 
short distance expansions) of the CG renormalization. These assumptions are 
consistent and amount to discard less relevant operators. These irrelevant op- 
erators can be identified within the method using replica of Section 3 where 
the above assumptions appear as standard in the RG of the m-component 
vector CG. For instance, the separation of the disorder into the two compo- 
nents V^-* and Vr corresponds in the replica method to the natural splitting 
in (11) between the vector fugacity local operator Y[n] (originating from Vr) 
and the off-diagonal replica Coulomb interaction Kaj^b (from V^). This will be 
further apparent on the equivalent Sine Gordon representation of the problem 
presented in Section 7. Accordingly, the definition of the independent local 
regions (and thus of the local disorder environments and of the detailed form 
of the distribution P{z~^,z~)) is clearly cutoff dependent. So is the detailed 
form of the fusion rule (26). However, universality of the physical results will 
be recovered in a remarkable way in Section 5, independently of the details 
of the cutoff procedure. As we will see, this is because, in the low tempera- 
ture limit, the above definitions and fusion rules capture correctly (to order 
P{z ~ 1)^) the rare events which dominates the physics. They correctly eval- 
uate the universal part of P{z~^,z~) (its tails in the low temperature region 
where they dominate the physics) while they also correctly describe the weak 
disorder regime at higher temperatures. 

Finally, we note that usual charge fusion between certain types of replica 
charges, represented on figure 4 corresponds, in the method without replica, 
to fusion of environments. 
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4-2 Expansion of physical quantities in the number of local regions ("number 
of points") 

To renormalize consistently the present model we need a method which can 
handle in a systematic way broad distributions of local fugacities. We have 
found such a method, which we now introduce. It is based on a systematic 
expansion of physical quantities in the number of independent points. It gen- 
eralises the conventional fugacity expansion in y of the pure case, but is more 
powerful. In effect, it amounts to a partial resummation of the conventional 
expansion. It is versatile since, as we will see, it yields back the conventional 
expansion in the pure case or at high temperature, but is also able to han- 
dle the broad distributions which arise at low temperature. The idea there is 
that only few rare regions (favorable to the charges) in each environment will 
dominate the observables and thus it becomes possible again to expand in the 
density of such rare regions. 

The idea underlying this expansion in the number of "independent points" is 
that the site-dependent fugacities associated with distinct "points" can be con- 
sidered as statistically independent. On a lattice, these "independent points" 
naturally correspond to the sites of the lattice, while in the continuum model, 
their definition is necessarily cut-off dependent. Upon coarse graining, we will 
be able to use this systematic expansion to renormalize consistently the fu- 
gacity distribution associated with each point. 

For definiteness, we will show how to construct the expansion in number of 
points on the free energy F[V] = In Z[V]. It can be performed in a 

given environment, keeping the full functional dependence in the set of po- 
tentials {Vr}. The construction can then be easily generalized to any physical 
quantity, such as arbitrary powers of the free energy (which yield an 

expansion of all moments averaging term by term over disorder) or all corre- 
lation functions of the field Ur which can be obtained from products of free 
energies . . . F"*:[T4] by differentiation with respect to the potentials K- 

In a second stage (next Section) we will use this expansion to justify the RG 
equation for the disorder distribution. 

We recall that the Coulomb gas model we consider is defined by its partition 
function 

Z[V]^1 + Y: E E /'^^^^r,-^Or,-r,-.+^J:.n^yr. (27) 

p>0 {ni,...np} rij^--j^rp 

Here and below, as in (7), all formulae can be extended to the continuum 
model by replacing discrete sums over distincts sites ri ^ r2 ■ ■ ■ ^ i"p by 
integrals with, e.g hard core conditions |ri — > and introducing the 
uniform fugacity y as was done in Section 2. Note that in the above expression 
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(28) wc do not make use of the decomposition (6) and K denotes the original 
disorder. 

The expansion in the number of points of the free energy has the form: 

F[V] = /(°) + E /S. [V] + E Al,r3 m + ... (28) 

ri^r2 ri^r2^r3 

where, in the disordered case, each term /rf-* depends on the values taken 
by the disorder potential V^. exactly and only at points ri, . . . r^. The explicit 
construction of these terms is given both for the disordered and pure case 
{y = 0) in Appendix G. Prom a practical point of view, the explicit expression 
of /if,^...rfc is in all cases: 

/S.r,=E(-l)'-' E Fr,^,-r.,[V] (29) 
/=0 h,...ii€[l,...k] 

where Fr^^^,,,ri^ [V] is the free energy of the Coulomb gas defined only on the set 
of I points r , . . . r (instead of the full lattice) and the summation is over all 
distinct subsets of the set ri . . . r^. The definition (28) is unambiguous, though 
subtle. Looking at the explicit expression (29) as a sum over smaller subset of 
points one could imagine adding other terms to the f^'^^ depending on less 
than k points. This is not possible in a global way, as the whole series must 
add up to the free energy, and the formula (29) enforces it order by order. We 
refer to the Appendix G for further details about the precise definition and 
construction. Note that the term /^^^ vanishes here because due to neutrality 
one cannot define a CG on a single site and that the expansion in the number 
of points of powers involves a rearrangment of the expansion of i^i^^]- 

It is important to stress that each term f^''^ of the above expansion corresponds 
to an infinite sum over terms of arbitrary high order in the conventional fugac- 
ity y (as defined in (9)). This can also be seen when setting disorder to zero, 
where we find that each f^'^''^ {k > 1) and (A; > 2) starts as y'^'^ plus an 

infinite number of additional higher order terms in y (see below). Indeed the 
expansion (28) corresponds to a complete resummation of the conventional 
fugacity expansion usually performed in Coulomb gas studies [9] except that 
it is usually performed on the partition function Z while here wc perform it 
on the free energy. The expansion (28) using the free energy is the appropriate 
expansion when the fugacities (or the core energies) are random and strongly 
site dependent with a broad distribution (e.g. most of the sites have ^ ~ 
except for a few which have z ~ 0(1)). Indeed, in this case the only small 
parameter is the probability P{z ~ 1) that a given point is favorable for a 
charge. Thus the k-th term of the expansion (28) is of order P{z ~ 1)*^, since 
it is associated with k independent points. We can thus consider (28) as a 
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perturbative expansion in the small parameter P{z ~ 1), valid in the XY 
phase, which replaces the conventional expansion in y. 

To compute the f^H^y.^ we thus have to consider a Coulomb gas defined by the 
partition function (27) (in the continuum limit), restricted to the system of 
points ri, . . .r^. We will in addition restrict ourselves to charges rir = ±1,0, 
as higher charges, examined later will be less relevant. Let us closely examine 
the lowest order terms k = 2 and k = 3. 

4-2.1 independent dipole approximation 

For k — 2we need only two points and the partition function (27) reads simply 

/I _ I \ "^-^-^ 

Zr„r,^l + W^r„r. = l+ (y "^^ + y^e^- "^^ ) (30) 

\ (^o J 

The three terms of the partition function (30) corresponds respectively to no 
charges or a dipole in {ri,r2}, and the two possible positions for the dipole. 
This results in the Boltzmann weigth for a dipole H^nrg- The first terms of 
(28) thus reads 

Al = -r'ln(l + W^r,rJ (31) 

Restriction of the expansion (28) to its first nonvanishing term, i.e F — 
X]ri,r2 /r?r2 corrcspouds to the so-called independent dipole approximation. 
This is the approximation on which the analysis of [23] was based. This ap- 
proximation, which neglects all interactions between dipoles, may seem at first 
sight to be a good enough approximation in the XY phase and at the transi- 
tion. Although this is reasonable in the pure XY model, this turns out to be 
incorrect here: by discarding the next term /^f^r2 rs one throws away crucial 
statistical correlations. Indeed, when renormalizing the distribution of local 
fugacities, we have to take into account correlations between dipoles induced 
by the disorder, which arise as follows. Suppose that the site ri is favorable to 
the creation of a -|- defect, i.e if 2;+ ~ 1 and r2 and are both favorable to 
creation of a — defect (while other neigbouring sites are unfavorable). Within 
the independent dipole approximation the dominant configuration would be 
to put both a dipole in ri,r2 and one on ri,r3 to take advantage of these 
three favorable sites. These are the configurations (4-6) in Fig. 6. However 
this configuration is forbidden because of the hard core constraint (and we 
have restricted to ±1 charges, higher ones being less favorable energetically). 
Thus we need to take into account the effective correlations between dipoles 
which arise because of rare favorable sites. This is done by considering the 
second term of the expansion /rf,r2,r3- Furthermore, as we will see below, con- 
sistent one loop renormalization requires to examine all terms in the expansion 
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(1) 



(2) 



(3) 




Fig. 6. Possible configurations for three sites: (1-3) are the only configurations 
physically allowed, (4-6) are present in the independent dipole approximation but 
forbidden in real model (see text). The bars denote pairs of ±1 charges 

(28) and how they change under coarse graining, and thus to go well beyond 
the independent dipole approximation. 

4-2.2 third order term and fusion 

Let us derive the explicit formula for f^f^^i rs- The partition function with three 
sites reads: 

F[Vr„Vr„ K3] = -P~' In (1 + Wr,r, + W^.rs + W^.r,) (32) 

However, since all terms of the expansion of /r^r2r3 in terms of W depend 
exactly and only on rir2r3, we have to substract to this free energy the terms 
depending on less than three sites, which can be identified as /if^2"'~/rir3+/r2r3- 
The final expression for the second term of the expansion of the free energy is 
thus 



&rs = In (1 + Wr,r, + W^.^s + W^^^s) (33) 

+P-^ [ln(l + Wr,r J + ln(l + Wr.rJ + ln(l + Wr.rs)] 

^ 1,(1 + W.,r,){l + Wr,rs){l + Wr^r^) J 

Let us first notice that now the three sites component of the expansion (28 ) re- 
stricted to the first two terms is exactly —P~^ In (1 -|- VF^ra + W^rirg + W^arg)- 
Thus adding this second term has cured the problem coming from the config- 
urations (4,5,6) of figure (6) which, as discussed above are not allowed. 

It is interesting to note that the term f^^l^^a present (and important) in the 

expansion of In Z, even though there are no actual configuration with three 
charges in a given environment (from neutrality). Indeed each term f'^''^ in 
the expansion (28) contains contributions from every even number (less or 
equal to k) of charges. As can be seen from above (fig 6) it takes into account 
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effective correlations between the distributions of fugacitics of three differents 
sites rir2r3 induced by the hard core constraints. By contrast, there would 
be no such term involving three sites (because of neutrality) in the similar 
expansion carried out on Z (conventional fugacity expansion). It was thus 
missed in previous studies, while in fact we show in the next Section that it 
gives rise to a crucial contribution to the renormalization of the distribution 
of disordered fugacities : the fusion of environments. 

The above defined expansion in number of independent points (a type of 
cluster-virial expansion) which we have illustrated on the first few terms, can 

be performed systematically to all orders. By coarse graining all terms of 
the expansion self-consistently, we will now obtain the renormalisation of the 
disorder distribution. 

4-3 Direct renormalisation on the free energy 

We now propose a renormalisation scheme based on the above expansion 
(28). From (4.1) one captures the relevant physics by defining random lo- 
cal charge fugacities at each point and splitting the disorder distribution 
as P[V] = Hr -P(^r^) ^ (6)- We know from the results of Section 

3 that one can renormalize the model by defining only three scale dependent 
quantities: the full local disorder distribution Pi{z'^,z~), the stiffness Ji and 
the second moment cr^Jf of the long range disorder V^. We now separately 
obtain the renormalization of each of these quantities, in a systematic pertur- 
bative expansion in the small parameter P{z ~ 1) deduced from the expansion 
in the number of points. There are three types of contributions as follows. 

4-3.1 Rescaling 

As can be seen from (7, 6) there is an explicit dependence in the cutoff oq in 
the expression of the interaction energy between charges and of the correla- 
tor of the long range part of the disorder . This dependence will appear 
in each term f^''^ of the expansion (28). Upon changing the cutoff this will 
results in contributions of order 0{dl) which can be absorbed by appropriate 
redefinitions as follows. 

The interaction term changes as: 

PJ ^ rir In ( ^— J ^pJj^nr In ( ^^^^ ) rir' - (5Jdl ^ (34) 

r^r' V / t^t' \ / r 

where we have used the neutrality condition = 0. The first term is the 
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interaction term with the rescaled cutoff while the additional term produces 
an additive contribution to the ±1 charge fugacity 

Similarly from (6) the correlator of can be rewritten as 

{V>{r) - V>{v')f^AaJ\hi ^— ^)(1 - ^^^^^r - r')) (35) 
+AaJ'^dl (1 - S^^°\r - r')) 

explicitly as the sum of a new long range disorder correlator with cutoff do — 
aoe^ and a short range disorder correlator (we have discarded terms of order 
0{dP)). Thus the original problem with cutoff Oq can be rewritten as one 
with cutoff So with (i) a new gaussian long range disorder with identical form 
of the correlator (6) with replaced by Sq (ii) a new local (short range) 
disorder v{r) v (r) + dv{r) with dv{r)dv{r') = 2aJ^dld^''°\r - r') since it is 
clear from (35) that when oq — > aoe^ the LR disorder produces an additive 
gaussian contribution dv to the SR disorder. Adding the two contributions we 
find that the change of cutoff produces the total rescaling contribution (25). 



4-3.2 Fusion of environments 

Having introduced the expansion (28) one can now coarse grain its continuum 
version by increasing the cutoff Oq — > = aoe°" and integrating over the 
corresponding degrees of freedom. Upon this increase two points r,, from 
the k points integral (/*^'^^) will be fused if they satisfy oq < |rj — r^l < aos^K 
This will produce a contribution which is an integral at scale Sq depending 
only on A; — 1 independent points and thus, by definition of the f'^''\ produces 
a correction of order 0{dl) to the k — 1 terms /(^^^). All these corrections can 
be reabsorbed into a correction to the fugacity distribution (together with an 
additive change to f^^^ the free energy contribution of all degrees of freedom 
which have been eliminated in the change of cutoff). 

We now illustrate how this works on the case k = 3, and indicate in the 
Appendix H how it works for arbitrary k. To lowest order in dl, this correction 
is independent of k and can be easily performed by considering the three points 
integral, using (33) with e.g Qq < |ri — < do — a^e^^ : 



Jao<|ri— r2|<ooe'*' ' ' Jao<|ri— r2|<aoe'* 



1 + 



r2 



- ln(l + PFr^rs) - ln(l + PTr^ra; 



(36) 



where we assume |ri — rsl > Oq, |r2 — ral > Oq. Upon coarse graining, the two 
points ri,r2 are fused to a single point f = (ri + r2)/2 and one obtains a 
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correction to /- j.^. 

Using tlic decomposition (6) of tiie disorder potential into a correlated com- 
ponent and a local part : Vr — V/' + Vr, we can rewrite the Boltzmann weigth 
for a dipole VFnra , defined in (30), as : 



/I _ I \ ' 

where we have used the definition — ye"^^^" for the local charge fugacities. 

Several simplifications now occur in evaluating (36). We first note that in 
(36) the integral is of order dl. Thus to find the leading correction of f~ll to 
order 0{dl), we only have to consider the integrand expanded in order in dl. 
To this order, from the correlator (6) of the disorder potential and using 
li"! — r2|/ao = 1 + 0{dl), we find that 

V> = V> + 0{dl) = T4> + 0{dl) (38) 

Thus the weight of the fused pair can be simplified as VFrir2 = (^^-^i^ +^Ii2-^^l)■ 
Similarly, VTrirs and Wr^r^ simplify to order in dl using (38) : 



W^r,r3 - f {zXz;^e^^^r-Vr%) + 4^-e^(^'^^^)) (39a) 

W^^ra - e^(^^-^^P + e^(^'^^^)) (39b) 

Using these simplifications, we can now rewrite the first term of (36): we obtain 
for the corresponding dipole weigth : 



-2/3J" 

1/1/ —1— 1/1/ / I I- fr. I \ 

W-L, 



r,rz 



1 + 



^ ( i^ti + ) ^rz ^f3(V> -V> ) (^n + ^ra)^!^ 

V 1 + ^r+i Z- + Z- Z+ 1 + Z+ Z- + Z- Z^^ 




It is now simple to verify that the change in the disorder averaged < Z^n^ra fr%2 [^1 
originating from the f^^^ term is correctly accounted for by following rules for 
the random fugacity variables. The last two terms of (36) produce, with the 
simplification (39a,39b), the rules 
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Prom (41), we find the rule corresponding to the first term of (36) : 

(Z^ Z±) , ^+ - (^r^l + ^rt) - _ (^n +^r2) /.q^ 

± T ^ri^r2 ' ^ri^r2 ^ ' ^ri^r2 ' ^ri^r2 

This is true for any point rs different from ri and r2. One can thus rewrite 
these fusion corrections as a correction to the (unnormalised) 'distribution' of 
local fugacities ^{z+,Z-) : 



di^z+,z^) = 02/ ^ ^z+,zj^z+,zj 



(43) 



X 



1 + z^^ Zj.^ + ^^j.-^ z^^ ^ 



The coefficient C2 comes from the integration over the relative position ri — r2. 
Thus from (43) and (25), we recover exactly the RG equation (18) for the 
function $. 

One can similarly check that the above rule correctly account for the correc- 
tions to k — 1 term from the k term upon change of cutoff. To be exhaustive, 
we must also consider a constant term in the expansion (28) of the free energy: 
/o(/) which satisfies difo ~ (ln(l + z'j^z'!_ + z'_z'\_)))i^^. This term corresponds 
to the free energy sum of all degrees of freedom which have been eliminated 
up to scale /. Indeed, P{zj^,Z-) (or $((^+,^_)) contains an average over all 
disorder configurations at smaller scales (all environments which have been 
eliminated). A full and systematic proof can in principle be made by consider- 
ing all averages of powers F^[V^] expanded in number of points and all fusions 
to order k giving corrections to order k — 1. We will not attempt it here, as the 
systematic study to all orders is much easy to perform within the (equivalent) 
replica formalism. 



4.3.3 Screening 

To derive the scaling behaviour of both the stiffness J and the strength of the 
correlated disorder , we consider the screening of the interaction and the 
correlation of the disorder between two infinitesimal test charges ei and 62 in 
the sample. To implement the study of this screening within our expansion 
in number of indcpendentpoints, we first define F[K 61,62] as the free energy 
of the disordered Coulomb gas defined in (7) with the two additional test 
charges ei, 62 fixed in ri, r2. These test charges interact with the other integer 
charges of the Coulomb gas, which screen the interaction between them, and 
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the correlator of the disorder V^. Prom F[V, ei, 62], one can define the screened 
interaction and disorder by 



/3G'^(ri,r2) 



1 d^FjV, 6^,62] 

2 deide2 

dF[V,ei,e2] 



61=62=0 



(44b) 



(44a) 



dei 



61=62=0 



At large distance we expect from renormalisability PG^ ~ /3 In ( | ri — r2 1 /oq) , 
which imply the definitions for the renormalized coupling constant and disor- 
der strength 



where the Fourier transform of the 2D laplacian G^^ has been defined after 
equation (2). These definitions can be transformed exactly using (7) into the 
relation 



Up to now, these are only standard definitions of the renormalised stiffness 
and disorder strength a^. From this one deduces the RG equations for J 
and a since by coarse graining (46a,46b) we obtain [25,24,23] 

TdiJ~^ = -2aoTT'^{nonR=a)c ; dia = -2aQTT^{no){nR=a) 

where the factor arises because we are dealing with correlations of the 
charge density. 

The novel and most tricky part is how to evaluate the right hand side in a 
systematic way. Our method is to expand these correlation functions in the 
number of independent points. This can be done in a systematic way using 
(28) and the definitions 




(45a) 
(45b) 



2a«4 = limGkVk«n 




(46b) 



(46a) 
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, ., , dF[V]dF\V] 

Note that while the expansion of the first term (connected correlation) in the 
number of points is readily obtained from the expansion (28), the expansion 
of the second term (disconnected correlation) is more subtle since we need to 
square the derivatives of (28) and rearrange it again in an expansion in the 
number of points. Fortunately, to the order we are working, we need only the 
two point term in the expansion. Thus we can restrict ourselves to terms with 
two points on both expansions. The expansion of the first term (47a) follows 
straightforwardly from (28) : 



dV^dVn r^Tr, dVodV^ ^^^^^^^ dV^dV^ 
and thus from (31) one finds: 



-13-' gWo,R 13-^ dW^^^dWo,^ 

W^CRWVk (1 + W^o,r)' dVo 5Vr 



since we can keep only terms depending of the two points (0 and R) in the 
last equality (the other contributions vanish). Similarly one finds that: 



dVo dVn {l + Wo,ny dVo dVji 



where the dipole weights Wo,n. have been defined in (30) 

These expressions further simplify since it is sufficient to work to lowest order 
in dl, i.e to order in dl. We can use that = Vq^ + 0{dl) since we need only 
to evaluate the correlation at |R| = Oo- This yields that VFo,r = -^o -^r + -^o -^r- 
One can also replace the derivatives with respect to the full disorder dy^ by 
derivatives with respect to the local disorder d^^ (see (6)) and thus we find: 
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al{no){nR=a) 



c 




(49b) 



(49a) 



(1 + Z^Zj^ + Z^Zq )2 



By using this expansion into (47a,47b) we find the screening RG equations 



Thus the expansion in the number of points has allowed us to derive consis- 
tent RG equations without using replicas. Note the difference with the pre- 
vious unsuccesful attempt in [23] to derive RG equations for the stiffness 
and disorder without rephcas. The authors of [23] were working in the in- 
dependent dipole approximation. Thus they computed the correlation func- 
tions as in the first line of (49) except that they kept the full expression 

(X]ri7^r2 dVo^ ) (I]r3^r4 Qy^'^ ) • This led to intractable expressions involving up 
to four independent points. The beauty of the present method is that we 
know that we can keep, in a consistent way to this order, only the terms in 
this expression which involve two points. 

4-4 Conclusion and connection with the replica method 

We have thus developed two RG procedures, one using replicas based on an 
expansion in the vector fugacity Y [n] which may have seemed somewhat for- 
mal, the other one, direct without replicas, as an expansion in the number 
of favorable regions P{z ~ 1) (at low temperature, while at higher tempera- 
ture it smoothly crosses over into the usual expansion in the uniform fugacity 
y =< z >). In fact these two approaches are equivalent. Indeed, we have 
checked using some combinatorics detailed in Appendix E, that the fugacity 
expansion in Y[n\ (13, 11) is identical term by term, in the limit m ^ to 
the expansion in number of points of the free energy (28). Thus, working with 
one or the other is equivalent It also shows that the limit of diluted vector 
charges physically corresponds to the limit of rare favorable regions. 



5 Analysis of the RG equations 

In this Section we analyse the RG equations derived in the previous Sections 
and we obtain the phase diagram and the critical behaviour of the XY model 

^ the RG procedure used in the previous Section to derive the screening equations 
without rephcas was found to be consistent with the more direct RG of the rephcated 
CG explained in Appendix C. 



(22). 
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with random phases. 

5.1 general analysis of the full RG equations 

Our task now is to study the closed set of RG equations (21,22) for the scale de- 
pendent distribution of fugacities -Pi(2;+, z_), the long range disorder strength 
(Ti and the stiffness J;. Since we study below the zero temperature limit, we 
have chosen to keep T fixed and renormalize J (only the combination K — 13 J 
flows) . 

We started by studying the set (21,22) numerically. Although we will not 
reproduce here the details we describe the main results. They confirmed the 
overall physical picture and led us to introduce an approximation, described 
below, which allows for an analytical solution. 

We found that at low T and for a smaller than some critical value Uc an XY 
phase exists (as in Fig. 1), where both the stiffness J and a converge to finite 
non zero values at large I: 

Ji^ Jr>0 (Ti^aR (50) 

In general, at low T we found that, starting at scale I — from (9) with a small 
initial averaged fugacity (and small local disorder) the distribution ^^(2;+, Z-) 
becomes broad and develops power law tails in the variables z_ which 
quickly extend up to fugacities z± ~ 0{1). However, in the XY phase (50) 
we found that the typical z± goes to zero and that the concentration of rare 
favorable regions, after an initial increase, ends up decreasing towards zero at 
large Z. It is useful to define the single fugacity distribution: 

Pi{z) = j dz+Pi{z+,z) = I dz_Pi{z,z_) (51) 

which does not satisfy a closed RG equation. Then we find in the XY phase 
that Pi{z ~ 1) decreases at large / (or cquivalcntly J^^^* Pii^) where z* is 
some arbitrary threshold). Thus the probability that either 2;+ or z^ is of 
order 1 decreases. On the other hand, for cr > cTc we find that Pi{z ~ 1) 
ends up increasing at large scale and the whole distribution Pi{z) ends up 
drifting towards increasing z values. This corresponds to the disordered phase 
where Ji vanishes at large scale {Jr = 0). The most interesting fiow occurs 
in the critical regime near the transition. There, if one interprets the fraction 
of favorable regions Pi{z ~ 1) as the perturbative parameter, the structure 
of the flow near a — ac is reminiscent of the RG flow a la Kosterhz-Thouless 
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with a separatrix, in the plane (-Pi(l). o") between the XY and the disordered 
phase. Accordingly, exactly at the transition, i.e on the separatrix (which, to 
be accurate, is a critical manifold in the space of distributions and couplings 
Pi{z+, J, cr) the distribution -P(-2+, becomes very broad and develops a 
fixed shape, with Pi{z ~ 1) slowly decreasing to zero (which makes the critical 
regime perturbative in Pi{l))- A schematic representation of the distribution 
is given in Fig. 7. 

Let us close this paragraph by noting that at higher T, by contrast, we found 
that the distribution remains peaked and that the overall picture becomes 
more consistent with considering a uniform average fugacity (as in [9]). The 
XY transition then occurs when this uniform average fugacity ceases to flow 
to zero at large scale. 

To quantify these (mostly numerical) observations, we now turn to the single 
fugacity approximation. 




Fig. 7. Schematic representation of the scale dependent behavior of the distribution 
P(z+,z_) in the plane of core energy variables {u,v) and the regions discussed in 
the text). The stifness Ji produces a convection towards the left (negative u axis) 
while at the same time the long range disorder a produces a diffusive spreading 
of the probability weight along the (local disorder) v axis (in both directions). For 
C7 > fJc diffusion wins and the distribution spreads. For smaller a < ac the weight 
of the distribution remains confined within the left quadrant 2:± < 1. In general it 
develops a two dimensional front solution. 

5.2 Single fugacity RG equation 

We now argue that, given the structure of the RG flow observed numerically, 
we can, with no loss of accuracy in all the regimes of interest (i.e within and 
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near the boundaries of the XY phase) approximate the full RG equations (21) 
for Pi{z^, z^) by a simpler equation for the single fugacity distribution Pi{z), 
and similarly for (22). 

Let us first focus on the low T regime, where the distribution Pi{z^,z^) is 
broad and the physics is dominated by rare favorable regions. For such a broad 
distribution the parameter which allows to organize perturbation theory is: 

Pi{l) = Pi{z ~ 1) ~ Pi{z+ ~ 1, z_ ~ 0) = Pi{z+ ~ 0, z_ ~ 1) (52) 

where we distinguish symbolically the rare configurations 2;± ~ 1 from the 
typical ones where z± 0. The first observation is that z+ and z^ are even 
more rarely simultaneously ~ 1. We note symbolically: 

Piil, 1) = Piiz+ ~ 1, z_ ~ 1) ~ Piil)^ (53) 

Indeed such configurations, absent from the start, are generated from the 
fusion term in (21). To generate from the fusion rule z± = {z'^ + z'^)/{l + 
z'^z'^ + z'_z'^) a configuration where both z+ ~ Z- ~ 0(1) one clearly needs 
at least that either z'_^ ^ z'!_ 1 or z'_ ^ z'l ^ 1, which is of probabihty of 
order -Pi(l)^ (symbolically one can write (9;P;(1,1) ~ Pii^)"^ + ..)• Thus, and 
this can be further checked, the estimate (53) is valid (the diffusion term does 
not change it). One can now inspect the corrections to the stiffness J and 
to the long range disorder a in (22) which come mainly from configurations 
where either ~ ~ 1 or z^, ~ ~ 1. This leading contribution is 
thus of order Pi{z ^ 1)^. Other fugacity configurations, e.g. with both 
z'j^ z'_ 1 contribute in (22) but, from the above estimate (53), their 
probability is ~ -Pi(l, 1)P/(1) ~ Pi{^T thus subleading in an expansion 

in PiiX). To summarize, one can see schematically the RG equations (22) as a 
correction of leading order Pi{lY to Ji and ai and the RG equation (21) as a 
correction to Pi{l) of order Pi{l) (diffusion term) and Pi{Vf (fusion term). 

Guided by these observations, we now approximate the full equation (21) by a 
RG equation for a function of a single variable which should correctly describe 
the behaviour of P(z+,2;_) around 2;± ~ 1. This approximation amounts to 
neglect the denominators in the fusion rule of (21). Indeed, in the spirit of 
our previous estimates, this denominator contributes only when 2;;^ ~ 2;" ~ 1 
(or z'_ z'^ 1) and then yields a contribution of order -P;(l)^ but which 
corrects only P;(l,l) since this fusion produces both a 2;+ ~ ~ 1. This 
corresponds to a correction of order Piil^ to -Pi(l) (which we neglect, since 
the leading order of the correction to Pi{V) coming from the fusion term is 
P;(l)^). This can be seen e.g. since Pz(l, 1) itself in turns contributes to Pi{V) 
to order Pz(l)^ (through integration over one of the two variables), or by more 
detailed arguments, not reproduced here. Neglecting the denominators gives 
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the equation for the probabihty distribution 



diP{z+, z_)=OP- 2P{z+, z_) (54) 
+2 {5{z^ - « + zimz_ - {z'_ + z'^)))^,^„ 

This equation can now be exphcitly integrated over one variable, say z_ and 
yields a closed equation^ for the single fugacity distribution -P/(-2+) defined in 
(51): 



diPiz) 



(3J{1 + zd,) + (y{l3Jf{l + zd,f] P{z) - 2P{z) (55) 
+2 / dz'dz"6(z - iz' + z"))Piz')P{z") 

Jz',z">0 



We can now obtain the corresponding equations for the corrections to J and cr. 
In (22) one keeps only the configurations corresponding to either z;^ ~ 2;" ~ 1 
or ~ 2;" ~ 1 but discard the much less probable configurations where all 
four fugacities are ~ 1. The expressions then nicely factor out in terms of 
Pi(z) and one obtains: 



C2 Jz' z'^ (1 + Z+Z 



where the additional factor of 2 counts the equivalent integral with + and — 
exchanged. 

We now summarize the obtained closed set of RG equations in terms of the 
variable u defined as: 



z = e^" (57) 

of distribution P{u)du = P{z)dz. Physically the random variable u can be 
interpreted as the random local core energy Ur = —Ec{r). One has: 



Note similarly that (54) admits solutions of the form P/(-z+, Z-) = Pl{z-^-)Pl{z-). 
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diP{u) 



(58a) 




+2 / P{u')P{u")5 (u 



8ci / /3e-^("'+"") \ 
^\(l + e-/^("'+""))2/ 



ci \(i + e-/5(«'+-"))2 




(58b) 



(58c) 



and thus the corrections to J and a involve only the convolution P P. 

We have thus justified, in the low temperature regime, that these approximate 
RG equations should describe the tails of the fugacity distribution exactly to 
the order -Pi(l)^ in the expansion in Pi{l) to which we are working. Indeed, 
since we have studied the RG in the previous Sections only to order y[n]^ (in 
the replica formulation, corresponding to -Pz(l)^ in the rare events formulation) 
it probably does not make sense at this stage to try to be more accurate, [f]. 

Finally, let us note that although we have focused until now on the low T 
regime it is rather obvious that the approximation of discarding the denomi- 
nators will be even more valid at higher temperature since in the conventional 
fugacity expansion (in a uniform averaged fugacity) these terms corresponds 
to Oiy^) terms while we work to Oiy'^). So these new RG equations interpolate 
all limits correctly. 



5.3 RG equation atT = 

One can now easily see on the form (58) that the RG equations admit a well 
defined T — > limit. In the equation (58a), the fusion rule u = ^ ln(e^" +e^" ) 
now becomes a max rule u = max(M', u") and (58a) transforms as 



Similarly, the function of u', u" in (58b, 58c) become respectively a delta and 
a theta function when T — >^ 0, so that in that limit one has: 



^ Although we believe the arguments based on an expansion in Pi{l) are correct, 
only an exhaustive analysis of the two dimensional front solutions of (21), clearly a 
complex task which goes beyond this paper, could confirm the exactness of the ap- 
proximation and the validity of the arguments based on organizing the perturbation 
theory using a single small parameter Pi{l). 



diPiiu) = Udu + aiJfdl) Pi - 2Pi + APi{u) j Piiu") (59) 

^ ^ Ju>u" 



39 



di{J-^) = ^ / duPi{u)Pi{-u) (60a) 



c 



2 



di{a) = ^- I du'du"Pi{u')Pi{u") (60b) 

C2 Ju'+u">0 

Note that the last integral exactly evaluates the probability to find two local 
regions with a total negative core energy, energetically favorable for a dipole, 
in agreement with the physical picture valid at low T. 

We thus obtain a close set of equations (59,60) which describes the scaling 
behaviour of the system at zero temperature. The equation for the distribution 
(59) can be conveniently simplified using the parametrization 



Gi{x) = / Pi{u)du (61) 

Jx-Ei 

with El — Jq JiidV . The function Gi{x) then satisfies 

\diG = ^-^dlG + G(l - G) 

In the case where a and J are /-independent, this equation is known as the 
Kolmogorov- Fisher equation and we will recall some results on this equation 
in the next Section. Before turning to its study, let us notice that the screen- 
ing equations (60) can be rewritten using this new parametrization simply 
substituting Pi{u) = —dxG{u + Ei). 

5.4 Parametrization atT>0 

Although (58) at finite temperature T > could be in principle studied di- 
rectly, it is much more convenient to introduce the generalization of the T — 
parametrization (61) as 

Gi{x) = l- [ P,(z) exp(-ze-^(^-^')) (62) 

Jz>0 

/oo _ 
du Pi{u) exp(-e'^("-^+^'^) 
-00 

with El = /q J{l')dl' . In the limit /3 — > cxd, this function reduces to the previous 
one (61). Using this new function and variable the integral equation (58a) 
can again be transformed exactly into a simpler differential equation which, 
interestingly, remains exactly the same as in the T = case: 

\diG = ^-^dlG + G(l - G) (63) 
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Only the initial condition explicitly depends on temperature (see below). 

Thus solving the KPP equation (63) allows us in principle to obtain the scale 
dependent fugacity distribution Pi{z) at any T. However the relation between 
this function Gi{x) and the distribution Pi{z) becomes much more involved at 
T > 0. As a result, and contrarily to the T = case, the screening equations 
do not admit a simple expression in term of Gi (x) . 

For any fixed I one can reconstruct the integer moments (z") p^ of the distri- 
bution Pi{z) by simply expanding the generating function in powers of e~^^ 
since: 

Gi{x) = Y.^—^{z-)p^e-^^^^~''^^ 

n=l 

Taking back this expansion in the KPP equation (63) and identifying the 
coefficient of the exponentials e~"^^, we find the exact RG equations for the 
moments (z") of the distribution Pi{z) : 

di{z") = (2 - npj + a{npjf) (z") (64) 
+ j dz'dz" [{z' + z"Y - {z'Y - {z"Y] Pi{z')Pi{z") 

Starting from a reasonable initial distribution with finite moments, the mo- 
ments remain finite for finite but, as we now discuss, increase quickly as 
I — > -|-oo. Let us examine the scaling dimension of the moments, neglecting 
for now the bilinear (fusion) term. We find that for fixed T and J, each suc- 
cessive moment (2;") diverges with the scale / when the long range disorder a 
becomes larger than a critical disorder a > a^"'^ = (T/nJ) — 2{T /nlY (see 
Fig. 8) (these values will be slightly renormalized by the screening equations 
but the conclusion will be similar). Putting back the fusion term simply im- 
plies that the moments with n > 2 also diverge when the first one does, as 
indicated in the Fig. 8. 

The divergence of the first moment yi = {z) p^ which can be identified with the 
uniform fugacity y of Rubinstein et al. yields the (incorrect) reentrant phase 
diagram of [9] where the XY phase is destroyed above the line a^^\T/J). 
Indeed it is already clear that the uniform fugacity approximation cannot 
work since we now see that higher moments diverge for even smaller disorder 
strengths. Thus even within the XY phase this result for the moments show 
that the distribution Pi{z) rapidly becomes broader and broader. Atypical sites 
where z is large appear and dominate the behaviour of the higher moments. 
Thus it becomes meaningless to study the scale dependence of the integer 
moments, but rather we must now consider the whole probability distribution 
Pi{z) and in particular understand its tail. This can be achieved with the help 
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Fig. 8. Critical disorder a^"'^ {T / J) above which the successive moments (z") diverge 
as the scale I — +oo. 

of the known solutions of the KPP equation which we briefly review in the 
next Section, before coming back to describing the scaling behaviour of Pi{z). 

It is interesting to note at this stage that the above KPP equation also arises in 
the problem of the directed polymer with quenched disorder on the Cayley Tree 
(DPCT) [29]. There the variable / corresponds to the number of generations 
and Pi{u) to the distribution of free energy —u = —Tin 2;. Thus we have 
demonstrated in an explicit and non trivial way, i.e at the level of the RG 
equations, that there are close connections between the two problems, which 
both exhibit a similar freezing transition. There are also notable differences 
between the two problems. For instance, while the diffusion coefficient D is 
constant in the DPCT studied in [29], in the disordered CG, Di depends on 
the scale and, in a self consistent manner, on the solution of the KPP equation 
itself, via the equations (56) which describes the physics of screening, absent 
in the DPCT. As a result, additional phase transitions exist here as will be 
detailed in Section 6. 

5. 5 The Kolmogorov-Petrovskii-Piscounov equation: some known useful prop- 
erties 

We recall in this Section some known facts on the Kolmogorov-Petrovskii- 
Piscounov (KPP) equation (also known as the Kolmogorov-Fisher equation) 
which we will need in the following Sections. The equation reads, in a general 
form: 



where the diffusion coefficient is constant, the function f{G) satisfies /(O) = 
/(I) = 0, / positive between and 1 and /'(O) = 1, /'(G) < 1 between 
and 1. The usual case corresponds to f{G) = G — G^. This equation has 
been applied to a wide range of problems, from chemistry to hydrodynamic 
instabilities or to the propagation of the Meissner phase into the normal phase 



^diG = DdlG + f{G) 



(65) 
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in a superconductor [35,36]. It is the prototype of equations describing the 
diffusive invasion of an unstable state by a stable one. This can be seen by 
writing it as a Landau equation diG = —dF/dG whose free energy F takes by 
construction its local maximum in G = (the unstable state) and its minimum 
in G = 1 (stable state). One usually chooses an initial condition Gi=q{x) 
monotonously decreasing from Gi=q{x — >• — oo) = 1 to Gi=q{x +oo) — 0. 
For a large class of initial conditions, the solutions of the KPP equations are 
known to converge uniformly towards traveling waves solutions of the form: 



Gi{x) — >• h{x — mi) (66) 

with h{x — > +oo) = and h{x — > — oo) = 1. However the qiiestion of the 
determination, given an initial condition, of the asymptotic traveling wave h 
and its velocity c = lim;^+oo dimi has been largely debated for KPP equations 
or for similar more complex non linear equations, and is still of current interest. 
It is known as the front velocity selection problem. 

It can be ilhistrated as follows. A family of possible front solutions exists, 
parametrized by the velocity c and noted hc{x — mi), as can be seen by substi- 
tuting (66) in (65). Constraints exist for the velocity. Indeed, one can linearize 
the KPP equation in the region ahead of the front for large positive x — mi 
where h is very small. As discussed below it is in fact this region which deter- 
mines the velocity and is universal. In this region one has: 



Dh" + ^h' + h^O (67) 
and thus h is a superposition of two exponentials e~'*(^~'^') with (see figure 9) 



c(/x) = 2(/x-i + Dii) (68) 



with c > c* = 4v D. The large x behaviour is dominated by the smaller fi and 
thus correspond to the left branch of the curve c(/i) in figure 9 when c > c*. 
In the marginal case where ji = jic = ^/\fD and c = c*, the two eigenvalues 
are degenerate and the front thus has the asymptotic behaviour for x — > -|-oo: 



h{x) - {ax + 6)6"''=^ (69) 

Let us now give the known results for the selection of the asymptotic front 
among the family of possible he- 

(i) Velocity selection. 

Although for more complex equations one relies on stability analysis and a 
" marginal stability criterion" [37] , in the case of the KPP equation a rigorous 
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4^ 

Fig. 9. Velocity of the front solution of KPP equation, as a function of its expo- 
nential decay ahead of the front. 

result is available. A theorem due to Bramson [38] shows that the asymptotic 
traveling wave is determined by the behaviour at a; — +oo of the initial 
condition Gi=q{x) in the following manner. If Gi=q{x) decays fast enough, as 
Gi=o{x) ~ e~^^ for x — >• +00 with jj, > Hc = ^/VD (or faster) theorem 
B of [38] states that Gi{x) uniformly converges towards the traveling wave 
solution hc*{x — mi) of velocity c = lim^^+oo 9;"^; — c* — A\fD. If Gi=q{x) 
decays slower, as Gi=q{x) ~ e~'^^ with /i < /i,. = 1/a/D then Gi{x) uniformly 
converges towards hc{^){x — mi) of velocity c{ji) continuously depending on /i 
and given by (68). The asymptotic velocity is thus given by c = A\/1D for steep 
enough initial condition fx> fXc= and c = 2(D/i + /x^^) otherwise. 

Moreover, the leading corrections to the velocity are also given by the theorem 
of Bramson and are independent of the function f{G) in (65). The correspond- 
ing position of the front m; is given by [38-40] 

for n> D-^ (70a) 
for /i = L>-5 (70b) 
for n < (70c) 

Note that cases (70a) and (70b) differ only by the velocity corrections but not 
by the asymptotic front shape which is he* in both case. 

Let us emphasize again the remarkable universality which arises in this prob- 
lem. Clearly the detailed shape of the asymptotic front depends on the detailed 
form of the non linear term f{G) in the KPP equation. However, the selection 
itself, the selected velocity, its corrections, and the tail of the selected front 
function h are all independent of f{G). It is also natural physically that it is 
the region in which the front penetrates (region ahead of the front x +00) 
which determines the selection. This universality has been explored further 
[40,41] and it has even been shown that the next leading corrections to dimi 
are also universal. Finally the marginal case which interpolates between (70a) 




mi (aI- - \nl + 0{l) 

^Vd (4l-l\nl + 0{l) 



--2iD 11 + 11-^)1 
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and (70b) has been also explored[^. The asymptotic front shape is the same 
for all these cases, but not the detailed finite / shape as we now discuss. 

(ii) Shape of the front for finite I 

The problem has also been studied for finite but large / [39-41]. We start 
with the case where Gi=q{x) decays fast for x +00 as Gi=o{x) ~ e~'^^ with 
/i > /ic = 1/Vd, most relevant for the following Sections. Then one must 
distinguish two regions in the traveling wave solution as illustrated in figure 
12. The central region is the "bulk" or "interior front" for |x — m/| fixed and 
finite. There, the shape of the front corresponds to its asymptotic form Gi{x) ~ 
hc*{x—mi) and the center moves as 9;m(/) = c* — ^^^ + 0{l~^) withe* = 4V^. 
Far ahead of this "interior front", x — rrii ^ 1, Gi{x) still decays faster that in 
the asymptotic solution hc*{x — mi). Thus there exist an intermediate region 
(which we call the "intermediate front region") which matches between the 
interior region and the region at infinity. In this intermediate region which 
corresponds to x — m/ ~ \^Dl one must take into account diffusion, and one 
can solve the linearized KPP equation without assuming a front solution, but 
assuming a scaling form. It is found [40,39] that in this scaling "far front" 
region x — ~ ^/DI, the solution Gi{x) behaves as[^: 



Gi{x) ^ gi{x - mi) (71a) 
gi{x) = A ( ^ + etc + 0(r 5) ) e-'^^^e"^ (71b) 



D 

Note that later on we will need to distinguish a region even further ahead of 
the front for x — mi ~ /. 

Let us close this Section by returning to the question of the dependence of the 
RG equation of the random XY model in the cutoff procedure. We can show 
that the cutoff procedure will determine the function f{G). Let us consider 
for instance again the limit of zero disorder, discussed in Appendix F. Since 
in that case D = one should have: 



-diGiix) = f{G) (72) 

with in that case 1 — Gi{x) =< exp(ye~'^'^^~^'^) >Pi{y) and P characterizes the 
"disorder" associated with the choice of cutoff procedure. An interesting choice 

^ It was found in [39] that when Gi=q{x) ~ ^.a^-iicx ^\^\^ _2 < a < one has: 
mi = y/D{^l-'^\nl + 0(1)) 

Similar results where obtained for the marginal cases G/=o ~ x'^e"'^"^ with —2 < 
a < with a scaling function noted Gi^((x — m{l)) / {2DI)) in [39] (for a = —2 one 

gets the above result, but the form for general a is more complicated). 
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corresponds to Pi{y) = S{y — yi) where yi satisfies the Kosterhtz Thouless RG 
equation diyi = (2 — K)yi. It is easy to see that in that case Gi{x) satisfies 
(72) with f{G) = -(1 - G) ln(l - G). This is further discussed in [16] 



6 XY phase and critical behaviour at low temperature 

6. 1 Phase diagram from the scale dependent fugacity distribution 

We now use the solutions of the KPP equation discussed in the previous Sec- 
tion to determine the phase diagram of the XY model with random phase shift. 
Let us first look for the XY phase where we expect that Ji and Di = \aiJi 
reach limits, respectively Jr and — \(JrJ'^ at large I. Thus in the XY 
phase at large I the KPP equation with constant D — Dr can be used. Pre- 
cise behaviour near the phase transition away from the XY phase, as well as 
intermediate scale dependence inside the XY phase a priori requires taking 
into account the / dependence of Di which will be done in the following Sec- 
tions. Note that I dependence of Ji itself results only in a shift that can always 
be trivially taken into account. 

First we note that the phase diagram will be entirely determined by the veloc- 
ity selected in the KPP equation. Indeed, we know from the previous section 
that Gi{x) converges to a traveling front solution Gi{x) — > hc{x — mi) of ve- 
locity c. The parametrization (62) then implies that the distribution Pi{z) 
of vortex fugacity for the random XY model, itself converges to a traveling 
front solution, more conveniently expressed in the random core energy variable 
u — \nz as: 

Pi{u) ^i^+oc P{u - Xi) , Xi=mi- El (73) 

where mi is given in (70) and Ei — Jq Ji/dl'. The center of the front of Pi{u), 
located in Utyp ~ + 0(1) corresponds to the maximum of the distribution 
P{u) (as can be easily seen on (61)) and to the typical values of the random 
variable u. The front shape of Pi{u) is simply related to the KPP front solution 
he through J^p{u) e-x.-p{e^^'^~^^) — 1 — hc{x). The asymptotic velocity of the 
front of Pi{u) is thus: 

- diE^yp = diutyp = diXi = dimi - Ji ^i^+oo c- Jr (74) 

The total velocity in the u variable is thus the KPP velocity minus the stiffness. 
The former comes from the spread of the distribution due to disorder, while the 
latter from the effect of interactions. In previous Sections we have explained 
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Fig. 10. Schematic representation of the connection between the sign of the total 
velocity of the front v = diXi = c — J and the increase or decrease of Pi{l)- 

that the XY phase corresponds to a decrease of the density of favorable regions 
Pi{l) = Pi{z 1) (i.e M ~ 0) and to the absence of topological defects at large 
scale. In particular, for J and D to reach finite asymptotic values Jr and D^, 
it is necessary that Pi{z ~ 1) decreases fast enough. Pi{l) can be estimated 
crudelyp^ as p{—Xi). Thus, as shown in fig. 10, when the total velocity diXi 
is negative (XY phase), the front moves to the left (large z or large u) and 
thus the probability of events z ~ 1 decreases while if diXi is positive, Pi{l) 
increases asymptotically (disordered phase). The XY phase thus corresponds 
to the region where the total velocity is negative and thus to: 



Jr>0 



(75) 



in which case the whole probability distribution of the core energy E^yp = —u, 
its typical and average values drift to +oo at larger scale. The transition line 
between the disordered and the XY phase can be located, in the plane Jr, a^, 
by finding the line where this relative velocity diXi vanishes. 

The phase diagram can now be obtained by determining the velocity c as a 
function of T and a. The crucial observation is that by construction (as is the 
case in [29]) the initial condition Gi=o{x) decays for large x as: 



-I3x 



(76) 



thus, since < z >Pp is finite, it decays exponentially and one can apply Bram- 
son's results (70) detailed in the previous Sections with the identification: 



/i = /3 



(77) 



Thus in effect it is the temperature which selects the velocity c. We find that, 
as depicted in fig 11: 



this estimate turns out to be accurate only near the transition as discussed below, 
but this does not change the conclusions 
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(i) for p < Pg — He — l/V-Dfl, i.e at high enough temperature T > Tg 
jR^Jcrji/2, the velocity continuously depends on temperature: 



c = c(/?) = 2iDnP + r = T (^2 + ^ j (78) 
and thus, in that regime, the XY phase exists for: 



9,X, = t(2-:^ + ^)<0 (79) 

which is exactly the condition which would be obtained from the averaged 
fugacity (at least when expressed in the renormalized parameters) and leads 
to the transition line of Rubinstein et al. as we have discussed earlier (see 
eq.(64)). 

(ii) the velocity of the front freezes (5 — (5g and for (3 > Pg, i.e at low 
temperature T <Tg, where it becomes temperature independent: 



c = c* = 4/d;= JrV8^ (80) 
and the total velocity of the front Pi{u) is now 

diXi^JniVS^-l) (81) 
Thus we obtain that below this freezing temperature at: 




the transition between the XY phase and the disordered phases occurs at: 
ctr = (Jc = |. For aji > ^ the system is unstable to the proliferation of 
topological defects induced by disorder. 



6.2 study of the XY phase 



To describe the XY phase it is important to understand the various regions 
of the scale dependent fugacity distribution Pi{z). Indeed the fugacity distri- 
bution also gives the distribution of the charge correlation functions in the 
Coulomb gas from the relations (49a). 
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disordered phase 




Fig. 11. Phase diagram as a function of the renormahzed disorder strength gr 
and the temperature T/ Jr. Tg corresponds to the freezing temperature and the 
transition hne for T < Tg is strictly horizontal when drawn in the renormahzed 
parameters (and slightly curved when drawn in the bare parameters). 

6.2.1 shape of the fugacity distribution in the XY phase 



far tail region ~J 




Fig. 12. Structure of the scale dependent fugacity distribution Pi{z) in the XY 
phase at low T and the several regions discussed in the text. 

In the XY phase one can describe the large scale behaviour of the fugacity 
distribution by neglecting the slow scale dependence of a and J (adiabatic 
approximation). The behaviour of Pi{u) can be obtained by studying Gi{x) = 
Gi{x + El) since one has 1 - Gi{x) = /„ Pi{u) exp(e^("-^)) with Ei = /^j dl'Ji^. 
At T = this simplifies into Pi{u) — —G'i{u). It is also useful to note that the 
Laplace transform of Pi{z) is Pi{s) — 1 — Gi{—T\ns). 

At low T in the XY phase Pi{z) becomes very broad and one must distinguish 
for large but finite I several different regions represented in fig. 12. Using the 
results of the preceding Sections on the KPP equation we now describe these 
regions in details: 

interior front: the first region corresponds to the bulk of the probability distri- 
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bution centered around the typical value of the fugacity Ztyp ~ ' . From the 
results for the front velocity (81,79) one thus obtains that in the XY phase at 
high temperature T > Tg = the typical renormalized fugacity decays 

with the scale as: 



which in that case is also the scaling of the average fugacity (2;)p,. On the 
other hand, below the freezing temperature T < Tg it decays as: 



zjyp ~ e-f^J{i-V8^)i (84) 

The bulk of the distribution is thus well described by the asymptotic front of 
the KPP equation as Gi{x) ~ h{x — Xi) and Pi{u) ~ p{u — Xi). Its precise 
shape is of course non universal as it depends on the details of the definition of 
the fugacity at the level of the cutoff (e.g. the function f{G) in (65)). However, 
since the physically interesting region z > Ztyp corresponds to the region ahead 
of the KPP front (even though the total velocity f = c — J is negative, see fig. 
12) universal results about the KPP equation can be used. In particular the 
near tail (i.e the scaling region u — Xi fixed but large) is universal and can be 
obtained since h{x) ~ xe^^^'^". Through Laplace inversion one gets that the 
fugacity distribution in that near tail region behaves as a power law: 




Pi{z)dz^[^ In (85) 



which is valid for < T < Tg and in the regime z/zj^^ fixed (as / grows) but 
large. The freezing which occurs at Tg thus concerns typical regions which, 
deep in the XY phase have a very small fugacity zf^^. It corresponds to the 
temperature at which the first moment of the distribution on the r.h.s. of (85) 
becomes infinite. Thus for T < Tg the true average fugacity {z)p^ ^ zj^^ and 
becomes dominated by rare local environments corresponding to values of z 
outside of the bulk of Pi{z) (where (85) is invalid). 

intermediate front: the second region is the "intermediate front region" (see 
fig. 12) and corresponds to u — Xi ~ Vl. There we know from (71a) that 
Giix) gi{x - Xi) with gi{y) ~ -^exp(-^ - ^). This region wiU be 
crucial to describe the critical behaviour in the following Section. 



far tail region: finally the most important region to describe the XY phase is 
the "far tail region" far ahead of the front with u — Xi ~ / (see fig. 12). Indeed, 
to obtain the renormalization of J and a (via the screening equations (56)) we 
are interested in the rare events 2; ~ 1 of small probability -Pi(l), but which 
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dominate the average correlations. These events correspond to the region of 
fixed u and thus to m — X; ~ —Xi {J — c)l. Fortunately, this region is so 
far ahead of the front that Gi{x) can be obtained with excellent accuracy by 
solving the linearized KPP equation: 

^diG = DdlG + G (86) 
by straightforward integration from from / = to /. This leads to: 



e ^J^Gi=o{x') (87) 
V Jo -^''"^ 

We also notice that in this regime Pi{u) can also be obtained explicitly. Indeed, 
since the relation between Pi{u) and Gi{x) is a simple linear (/-dependent) 
convolution it is straighforward to show (e.g. via x Fourier space) that Gi{x) 
obeying the linearized KPP equation (86) is equivalent to Pi{u) satisfying: 




diP{u) = [Jdu + oJ^dt) P + 2P (88) 

This linearized equation was not entirely obvious to guess directly from (58) 
and justifies, even in this simplest regime, the detour through the rigorous 
results for the KPP equation. Using Ei = Jq dl'Ji/ it yields: 



^21 



Piu) 



Sn /o Di,dl 



f 8 r D,,di' fj / / 

= / e Jo I Pq{^u 

I Ju' 



^9) 



It is interesting to note that for fixed x and u, Gi{x) and Pi{u) can be estimated 
in the large / limit as simple exponentials: 



with C = j^i e^'^J Pi=q{x') and an identical expression for Gi{x = u) with a 

x' 

different prefactor C = J^, c^^Gi^q^x'). These forms, which are valid in the 
region of fixed fixed u ~ 0(1) (and fixed x ~ 0(1) respectively) will be useful 
to estimate the renormalization of J and a belowf^ 

they were obtained from (87) and (89) by approximating exp(— (x^ + 
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Thus at low temperature we find the following decay of the probability of rare 
favorable local environments: 



Pi{z ~ 1) ~ A/-i/2g2i(i-g^) (91) 



with A = C / y4jw\P. We also find that the distribution Pi{z) has an algebraic 
tail at low temperature in the region 2; ~ 1 (i.e. z fixed as I 00) as: 



P,{z)^m)-^ , T* = 2aJ (92) 

with is a different power law behaviouip^ than the one which characterizes 
typical fugacities (85). These two different power law behaviours are repre- 
sented in figure 12. 

Finally, one can check that the three regions match properly and thus we have 
a fairly complete description of Pi{u). For instance, using the expression valid 
in the region (ii) for fixed x at large I one gets: 



Gi{x) ^ giix - Xi) ^i^^ y^e^ ^^yij+^i) 

= , e 

which always gives the result (90) up to the logarithmic corrections / prefactors 
factors, and reproduces even the / prefactors correctly for Xi/l 0, which is 
when we expect the matching to become exactp^. 



6.2.2 Screening and correlations in the XY phase 

To study the screening equations (56) in the XY phase we need the distribution 
Pfiud) = Pi*uPi of dipole core energy Ud = u' + u" . Fortunately in this regime 
it can be computed simply from (89) and is given simply by: 

x''^)/(8Dl)) ~ 1, and similarly for u (to be precise the D and J which appear 
here are l/l Jq Di and 1/lJ^Ji respectively). They are valid only as long as the 
integrals which defines C and C are convergent. Since Gi=q{x') ~ e'^^^' for large 
positive x' the estimate for Gi{x') is valid only as long as T < T* = 2(tJ. for T > T* 
one must instead perform saddle point estimate, see below 

note that for > 7"* the averaged fugacity {z)p^ is controlled by z <C 1, while for 
T < r* is controlled by z » 1 

this type of matching has been used to derive the | coefficient in the logarithmic 
correction to the velocity [40]. 
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Pt{n)= I \ e ^^lo-^-" P,\u')^C"^===e-^ (93) 



the last equality being valid for fixed u and Z — > +00, and C" — /„, exp{-^)pQ{u') 
in terms of the initial dipole core energy distribution. Substituting this last 
form in (56) we obtain, in the large I limit: 



dJj-^) = ^ C^^ / = Cjmi-^ Pi(lf (94a) 

dia = ^ , C" / di^^- ^ = CslMl^ Piil? (94b) 

with Cs — 8ciC" J-v^ttct/ (-\/2c2C^) a constant roughly independent of the tem- 
perature. The two above integrals 



T 



TT 



l-^)7r 



T* sin(^7r) 



sm(^7rj 



(95) 



are convergent respectively only for T < T* and T < 2T*. in which cases it 
is indeed legitimate to replace Pf{u) by its above asymptotic form. Note that 
the T = limit is well defined since in that case Ij = I and I„ = T* = 2a J. 

At higher temperature T > T* the full front solution controls the renormal- 
ization of J. Using the simple above gaussian form yields that the u integral 
is dominated by the saddle point u — u' — 2Jl + 4aJ^l which gives: 



di{J-') ^ e^^''-*+^^ T>T* (96) 

which corresponds for T > T*, to the behaviour of < 2; >^ (and for T > Tg of 
zfj^p), as expected. A similar, though more involved, analysis can be performed 
for a. 

Thus we have obtained the equations (94) for the renormalization of J and a 
in the XY phase. Since we have shown that Pz(l) decreases exponentially as 
(91) we conclude that J and a reach their finite limits Jj? > and ur < +00 
as power laws of the systems size. 

This analysis also yields the full distribution of the correlation function of the 
charges in the XY phase. Indeed, let us recall that: 
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-aQ{nonR=ao)c 
-ao(no)(nR: 



(1 + ZoZ^ + 4^0 )^ 
2 



1^0 ^R 



Zt3 Z, 



R^O 



^1 + ^0 ^R + 



(97a) 
(97b) 



Thus by following the RG up to the scale I* — ln(i?/ao) = ln(ao/ao) we can 
obtain from P;« {z) the distribution of the charge correlations at large R. 

For instance, from the above and (84) one finds the following decay of the 
typical thermal and disorder correlations: 



-{no)(jiR) 



typ 



typ 



aoy+¥(i-^) 

'rj 

aoy+¥(i-^) 
'rj 



(98a) 
(98b) 



for (7 < CTc and T < Tg. For T > Tg one has instead {nonR)c\typ ^ (ao/-R)^'^/^~^''•^^/^^ 

The averaged moments can be obtained as above by substituting the exponen- 
tial form of the distribution Pf*{u) for m ~ 0. Performing the corresponding 
integrals, one gets: 



(99a) 
(99b) 



with 



BpiT) 



C" 



C" 



-.T 



-.T 



r 




r 


P- 


jl' 
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[2p] 






r 


" T ' 
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2p- 


_ X 





VSW^^ T[2p] 
These formulae are valid only at low enough temperature T <T* (p) . 



(100a) 
(100b) 



Finally, although we have not attempted a precise RG calculation of the XY 
order correlation functions, the following behaviour should hold in the quasi 
ordered XY phase: 



(101a) 
(101b) 
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with 77 — + j^) and rj = ^(Jr- At the zero temperature transition point 
cr — cTc — 1/8, T — 0, the value of these exponents are universal rj — rj — 1/16. 



6.3 Critical behaviour at zero temperature 



We now study the transition from the XY to the disordered phases at low 
temperature at ctr = 1/8. From the previous Sections we know that it occurs 
when the total velocity of the front of the distribution Pi{u) vanishes, i.e the 
critical region is defined by 



diXi = J{V8a - 1) ^ (102) 

in a large / regime. While in the XY phase the physics was dominated by the 
far tail of the traveling front -u — X/ ~ / in the critical regime it is the interior 
front u — Xi ~ 0(1), as well as the intermediate front region u — Xi \/l 
(see fig. 12) which controls the transition. Thus the correct description of the 
transition requires the knowledge of the KPP physics in an essential way, and 
is thus entirely novel. 

For simplicity, we will only present the analysis at T = 0. We will also work 
to first order in cr — Uc- We first assume that the coefficient D — \cfJ'^ varies 
sufficiently slowly near the transition so that the results from the KPP equa- 
tion with a constant D can be used. This assumption will be self consistently 
verified at the end. 

Near the transition the center of the front is located u — Xi with Xi fa 
(4-\/D — J)l — ^\/D\nl + Xq as indicated by (70). Thus in the critical regime 
one still has Xi —00, although logarithmically in / exactly at the transition. 
Thus this critical regime can still be studied pcrturbativcly, as Pi{l) remains 
very small. The front velocity has the following scale dependence: 



diXi = Ay/D -J- + h.o.t (103) 

where the h.o.t. contains the universal [40] 0{l~^^^) subdominant corrections 

to velocity in the KPP equation as well as additional subdominant corrections 
originating from the slowly varying Di. In the critical region at T = we can 
use the KPP front solution Pi{u) = pi{u — Xi) with, from (71a): 



pi{x) = -g'i{x) A^e sdi 



(104) 
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an expression valid as long as x/l — > 0. The RG equations for J and a in the 
critical region thus read: 



di{J-^) = ^ / dupi{u - Xi)pi(-u - Xi) (105a) 

C2 J 



dia^^f Pi{u)pi{u') (105b) 

C2 Ju+u'>-2Xi 



We need to evaluate these expressions for Xi large and negative (typically 
either as ~ — In Z on the critical manifold or as ~ (a — ac)l very close to it). At 
criticality they behave approximately as exp(2Xi/y/D), which can be guessed 
by setting -u ~ in the first expression. A more accurate estimate of the above 
integrals is performed in the Appendix J. The end result is that, in the regime 
of interest (where we can discard terms of order {a — acY) one has: 



a,(j-i) ^ _^Xfe2^'/^ , dia CXfe'""^/^ (106) 



where C is a constant. Introducing the small parameter: 



gi^e^ (107) 
the density of favorable regions reads: 



^ Xi x^ 

Pi{l) ~ A^eVhe-sii ^ XiQi (108) 

since we can discard terms of order {a — (7c) ^. From (103) one finds that gi 
satifies precisely the RG flow equation: 



dig = (l6{a - ae) - ^ + h.o.?j g (109) 

to lowest order in a — ac- Note the | universal factor which arises from the 
universal velocity corrections in the KPP equation (70). 

A natural choice of parameters to describe the universal behaviour around the 
transition is thus g (which up to logarithmic corrections is equal to the density 
of favorable sites Pi{l)) and a which satisfies: 

dia^c(\n^\ / (110) 
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These two equations (109, 110) form our complete set of RG equations pro- 
jected on the plane a,g. They are somewhat analogous to the one describing 
a Kosterlitz-Thouless type transition, with an important difference. Here one 
readily finds that the separatrix is vertical. Introducing the deviation from 
criticality, e = cTc — (Too (where aoo = c^r) one has for e = the flow: 



9i 



, a = a^-7(ln/fr2 (111) 



and thus satisfy a^o — cr^ ~ g^^^ (\n.{l / gi)Y . Thus our RG equations yields 
gi ~ at criticality|__] and a correlation length: 



e ~ exp(- -) (112) 

since starting away from criticality e > one findsQ3 gi ~ l-^e^^'K This 
critical behaviour correspond to a new universality class which is different 
from Kosterlitz-Thouless and from the prediction of [24,25]. Note the crucial 
role of the | universal factor. Replacing it by any number less than 1 would 
have led to usual KT behaviour. 



We can now check that the variations of D should be unimportant at the 
transition. Indeed one find^^ that D^o — Di ^ (In/)^/^^ at most. On the 
other hand we can estimate that if Doo — Di varies faster than dirrii ~ 1// the 
KPP results should not be affected (see the scaling with / of all terms in e.g. 
Eq (A7) in [39]) 

Let us close by noting again how universality appears in the derivation of the 
critical behaviour. Although most details of the fugacity distribution P{z±) 
(e.g. its bulk, the fusion rule..) depend on the cutoff procedure the univer- 
sality in the XY transition appears in a remarkable way. It arises from the 
independence [38-40] of the velocity, the velocity corrections, and front tail on 
the precise form f[G] of the non linear term in (65). This gives us confidence 
in the method developped here. 



this indicate that at T = {nonR) ~ R~^{lnR)~^ at criticality 
the critical finite size corrections to a read fioo — ai = e^F[d] with F[x] = 
J"*""" dxx~^e~^'^^ up to logarithmic corrections. 

^^note the following interesting cancellation in the variations of Di. Writing dia = 
Ai[Xi] as a function of Xi, we get that di{J^^) = ^A'[Xi] since Xi appears only 
explicitly in the bound of the integral in (105). Thus one finds diD = ^j'^A[Xi]{l — 
Ja-^lnA[Xi]). Using that using that -^In A[Xi] = 2/\/I) + o(l) one gets that 
the leading variations cancel at the transition point a = 1/8. 
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7 Equivalent Sine Gordon model 



As is well known the Coulomb gas can also be equivalently formulated as 
a Sine Gordon model (see e.g [42]). In this section we identify the random 
version of the Sine Gordon model to which our analysis applies. 

The scalar Sine Gordon model, of partition function Zgg = J D(f)e~^^'^^ is 
defined, in the absence of disorder by the action: 



S{4>) = I (fv{^{WA f-gcosm (113) 

where 0(r) G [— cxd, +oo\. As is well known it is equivalent to a Coulomb gas 
since, expanding in g one has: 

+ 0O 1 

^^. = E7^/"((/^'rcos(0(r)))2Oo (114) 



p=0 



hoo 

^i{<t>ira)-<t>ira))\^^ 



P = P- Q = l 



where (..)o denotes averages with respect only to the quadratic gradient part. 
The interaction G{r) = r(0) — r(r) ~ In |r| has been defined in (2). The above 
partition sums involve only p neutral pairs of dipoles since, in the large size 
limit r(0) —>■ oo. Finally 7 = 1, ..2p in the last line involves a summation over 
all distinct charges with = 1 for 7 = 1, ..p and — —1 for 7 = p + 1, ..2p. 

Let us turn to the disordered version of the model. To reproduce the bare 
version (7) of the model, one must first add a short range correlated random 
imaginary field as follows: 



5(0, A) = J o^2r(^_i_(Vr0)2-^i-A•Vr0-^cos(0)) (115) 

with correlations AgA_g = na. Since it imposes now that (iVr0)o = 2i^A, 
each factor exp(mr0(r)) in the g expansion in (114) yields an additional 
exp{nrPV{r)) where P'^VgV^q — and thus reproduces the bare CG ver- 

sion (7) of our model. 

Note that the above model still contains a uniform "fugacity" g. It thus cor- 
responds to the version studied in [9] by Rubinstein et al. However we know 
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that this cannot be the correct form under renormahzation. The first obvious 
idea is to generahze g — ^ ^'(r), i-c a disordered fugacity (the above expansion 
can be immediately generahzed to this case). This term wiU be generated, but 
is not the end of the story. Indeed let us consider the symmetries of the above 
action 5'(0) (even in the presence of a g{v)). When A(r) = the action is 
real and invariant through (^(r) — 0(r). In the presence of the random field 
this symmetry is broken (as 0(r) acquires a non zero, disorder environment 
dependent average) and the action is complex. Thus nothing prevents that 
under coarse graining the action will become: 

/ (s^^^'-'^^' ~ ■ ^^'^ ~ ^+i^y^''^ - ^-(r)e-^^«) (116) 

and this is indeed precisely what we have found in the CG formulation, i.e 
each sign of charge acquires a different local random fugacity. 

There is however a symmetry constraint on the distribution of the local ran- 
dom fugacities. Indeed in the above bare model (115) the full partition sum 
is real, as there is still a statistical symmetry: 



5(0, A) =5(0, -A)* (117) 

for any configuration 0(r) and environment A(r). Since the probabilities of 
A(r) and — A(r) are the same, all physical averages will be real. In the coarse 
grained model (116), since we have that 



5(0, A, z+, z_) = 5(0, -A, z^y (118) 

the probability (over environments) of the random fugacity disorder configu- 
ration {2;_|_(r), ^-(r)} should be equal to the probability of {^_(r), Zj^{y)}. 

We can now check that this random Sine Gordon model (116) is indeed equiva- 
lent to the coarse grained disordered CG considered in this paper by expanding 
its the partition sum in a given environment, which yields: 



+ 00 \ P f 

Zsg{A,z+,Z-) = X! n / d^rad^iaZ+{ra)z-{Yo) 

p=0 P- a=l 



X exp K^^n^n^G{v^ - r^) /3l^^(r^)n^ (119) 

and thus the A random field disorder is associated to the long range part 
of the disorder in the CG (see (6)) while the random couplings constants z±{y) 
in (116) are associated to the local random fugacities in the CG. 
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Similarly one can establish a correspondence directly on the replicated versions 
of both models. Replicating the above SG model and averaging over disorder 
yields (the limit m — > is implicit) : 



■' a 

X /e/<^'rz+(r)^^e^^a(-)-z_(r)^^e-^^<^(-)\ .^20) 

\ / Z+,Z- 

Expandcd to first order and treating the -2+(r) and -2-(r) as uncorrelated in 
space and reexponentiating yields the replicated SG model defined as: 

^ = / n ^'^"^-^""^'^"^ (121a) 
SrepM = / dh I ^ ^(l5„, + a)V<j>a ' V (t>A - ^ y[n]e-"^ (121b) 

\ a,h ) n^O 

where n - = n°-(^a is the scalar product in replica space and y[n] are anal- 
ogous to the vector fugacities introduced previously in the replicated vector 
CG model. Although the bare model obtained from (120) contains only single 
component replicated charges, all multicomponent charges will be generated 
upon coarse graining, as in the rephcated vector CG. As we have seen all these 
charges should be taken into account and thus the generic replicated SG model 
should contain all possible y[n] with n 7^ 0. 

This SG model can also be studied using RG, either in its replicated form 
(121b), or directly (116), very similarly to the vector CG studied in this pa- 
per. Although variations in definitions of the fugacities y[n] and different 
cutoff procedures can induce some irrelevant differences in the details of the 
renormalization (and different RG equations), the two models have the same 
physics. The Sine gordon formulation has several advantages, such as exhibit- 
ing by construction the decomposition of the disorder in two physically very 
different components (see 6). It is also more amenable to replica variational 
methods than the CG, left for future study. 



8 comparison with previous approaches 



Let us compare our method and results with the work of Scheidl [24]. In 
this work the multicomponent charges (restricted to 0, ±1 in each replica) 
were considered, but the fusion was not taken into account. Also only dipole 
fugacities were introduced. 
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We can recover the RG equations and results of Scheidl within our approach 
by (i) artificially setting the fusion coefficient to in (18) (ii) assuming a 
log- normal distribution P = / $ for the fugacities (which is consistent only 
when fusion is neglected) (iii) defining random "dipole fugacities" as z\ = 
z'_^_z'!_ and = z'_z'l (or equivalently u = u' + u") identical to the "dummy 
gaussian variables" introduced by Scheidl. Within our diffusion formalism (see 
Section(3.3)) we find that the norm / $ diverges exponentially, which allows to 
recovers the extra factors e^' appearing in Scheidl's screening equations. Going 
from our diffusion formalism back to the replica formulation yields back the 
RG equations of [24]. 

In presence of a fusion term C2 > 0, the equations become a priori very dif- 
ferent. The fugacity distribution does not remain log-normal as we do not 
assume a priori the form of the distribution. Interestingly, although the log- 
normal does not reproduce correctly the true distribution of fugacities (and 
misses connections such as the one with the freezing of the DPCT via the KPP 
equation), some of our results deep in the XY phase, such as the renormaliza- 
tion of K and a (see (94)), agree with the one of Scheidl. It was not obvious 
a priori that the approach without fusion in the XY phase did not miss extra 
relevant physics in this model, and indeed near the transition fusion appears 
to be crucial and must be taken into account. 



9 Conclusion 

To conclude, we have constructed in this paper a novel renormalization group 
method which allows to study perturbatively, and in a consistent way, a 
large class of disordered models which can be formulated as two dimensional 
Coulomb gases with quenched disorder. We have applied it specifically to the 
XY model with quenched random phases. We have obtained the phase dia- 
gram for this model, confirmed the existence of a low temperature XY phase, 
and elucidated the critical behaviour at the transition where topological de- 
fects proliferate. It would be interesting to check our predictions in numerical 
simulations!^ 

The present RG method is not based on the conventional perturbative ex- 
pansion in a vortex fugacity y spatially uniform over the system, which, as 
we have shown, is only justified for pure models or for disordered models at 
high enough temperature. Instead, it is constructed by first defining the local 
random vortex fugacity (or core energy) and then following its full proba- 
bility distribution under coarse graining. Below a freezing temperature this 

in particular the full distribution of local fugacities for finite sizes should be 
directly measurable in simulations. 
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distribution becomes very broad and cannot be followed by conventional CG 
methods. Our renormalization procedure allows to follow this broad distribu- 
tion in a controlled way, by defining the new perturbative parameter as the 
concentration P;(2; ~ 1) of rare regions favorable to vortices. This running pa- 
rameter flows to zero in the XY phase and at the transition (marginal flow). 
Hence both can be studied perturbatively. The underlying physical picture 
obtained here is that the transition is controlled, as the scale increases, by 
the proliferation of vortices in less and less rare favorable regions. We find 
that it has some features reminiscent of a Kosterlitz-Thouless transition, with 
important differences, such as the scaling of the correlation length. 

To derive the RG equation for the distribution of vortex fugacities, we have 
introduced two equivalent methods. One is based on the replicated vector 
Coulomb gas version of the model, and in an expansion in the vector fugac- 
ities. The second one is direct, with no use of replicas, and is based on a 
systematic expansion of all physical quantities in the number of points, i.c in 
independent local regions and thus, in the end, in powers of the concentration 
Pi{z ~ 1) of rare favorable regions. As we have shown these two methods 
are fully equivalent: the first one being more systematic and the second one 
allowing for a clear physical understanding of the problem in terms of proba- 
bility distributions. These two expansion methods are highly non perturbative 
in the original uniform fugacity variable y. Since they are constructed from 
charge fugacities they can be made fully consistent (by contrast with previous 
approaches based on dipole fugacities). 

Our method sheds light on the broader issue of universality in random systems. 
The spirit of the RG method is that at large scale most information about the 
system is irrelevant and can be discarded. In constructing our RG procedure 
we have first shown, using the very special properties of logarithmic interac- 
tions, that it is enough to follow, in addition to the two parameters K and a, 
the distribution of only one or two local (i.e uncorrelated) random variables 
P;(^+,2;_). At this stage it is clear that we still keep too much information. 
Indeed we have found that the precise form the non linear RG equation obeyed 
by Pi{z^, zJ) as well as the detailed shape of this distribution are largely cutoff 
dependent. However, this complicated looking RG equation can be generally 
recast, up to irrelevant terms, as a well known non linear front propagation of 
the KPP type. Using the known remarkable universality property of this type 
of non linear equations, we found that all the information needed to describe 
the universal properties of the XY phase and of the transition is indeed in- 
dependent of the detailed shape of the fugacity distribution, or of the precise 
form of its RG equation. Only its tails, and the finite size corrections to the 
front velocity seem to be needed to determine the physical quantities and the 
critical behaviour. 

Since we are following a full distribution of local disorder, the present method 
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could be termed a functional RG. We can indeed draw some parallel between 
this functional RG procedure for the disordered CG and two other known ex- 
amples where the universal behaviour of a disordered finite-range systemp^ is 
extracted from a functional RG equation. The first one is the asymptotically 
exact real space RG in d = 1 [43,44] well suited to "infinite disorder" fixed 
points. As emphasized in [16] there are indeed similarities when treating the 
single vortex, d = 1 version of the present model. The method of [44] can be 
applied for the Sinai potential which has correlations growing as a power of 
distance. For the present case of weaker, logarithmic correlations, these meth- 
ods can also be applied in principle (at least around zero temperature where 
disorder is still very broad) but become very hard to work out analytically. The 
other example of functional RG appears in a dimensional expansion for the 
problem of an interface in a random potential [45] , which has infinite number 
of marginal directions at the upper critical dimension d = dc- An important 
question is whether in the problem studied in this paper there is also an infi- 
nite number of marginal directions. As discussed above, the results extracted 
from the KPP equation seem to suggest that a smaller amount of information 
than the exact full distribution may be needed here. Thus one can speculate 
that the critical theory studied here could be equivalently formulated as a 
more conventional field theory, yet to be identified, with a small number of 
marginal or relevant operators. In any case the RG method developed here 
should provide a physically transparent guide to study the system. Given the 
wide applications of Coulomb gases in two dimensions it is likely that other 
two dimensional disordered models can be studied using methods similar to 
the one introduced here. 

Finally, another outcome of the present work has been to unveil some interest- 
ing connections between the renormalization of a disordered system and the 
universal features of the propagation of invading fronts in non linear systems. 
The existence of intriguing relations between freezing transitions and velocity 
selection in non linear fronts was noticed previously by Derrida and Spohn 
[29] in their study of the DPCT. Here, we found an even deeper and puz- 
zling connection, betwen the "universality" in these selection mecanisms and 
the universality in the critical phenomena captured by the renormalization 
groupp'^. Since attempts have been made to construct renormalization meth- 
ods in order to extract the universal features of such non linear equations [46], 
this suggests that a common framework could be developed in connexion with 
two dimensional disordered models. 



^ and statistically translational invariant by opposition to e.g. Migdal Kadanoff RG, 
or Cayley tree recursion relations, which are exact only on very special hierachical 
lattices. 

'^^ the dimensions of operators, in the RG sense, being directly related to the selected 
velocity of non linear fronts 
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A Disordered lattice Coulomb Gas 



In this appendix we derive, by exact transformations, a lattice disordered 
Coulomb gas formulation of the Villain form of (1), extending to this disor- 
dered case the approach of Kadanoff [5]. An alternative route to [5] in the 
pure case in the method used in [32]. For simplicity wc firt turn to the Villain 
version of (1) before making the duality transformations, although the inverse 
procedure could be used (see [32]). This Villain model corresponding to (2) is 
[47]: 



Villain 



27r 



exp 



(A.l) 



where K = f3J and we used the notation h^x/y = U,i+e^/y (see fig. A.l). This 
Villain partition function indeed corresponds to a Z-gauge theory, since the 
action is invariant under 6i ^ 0[ = 9i + 2nqi, lij — > I'^j = kj + qi — qj with q^ e Z. 
Choosing the gauge field qi such that for each horizontal link l^^^ — (li+ex ~ Qi 
{i.e I'ix = 0) (free boundary conditions along x) we obtain a partition function 
of a gaussian field 0, = ^,/27r (Ez,. /X" f = IZo d(l>i): 



^Villain = 2^ H / "01^ ^ ' 



with the notation for the discrete gradient Vi^yCp = (pi+iy 
this gaussian field yields the Coulomb gas action 



j. Integrating over 



-^Villain — 6 



-PJ Ei J [('-% .y-^i,v)- 2^ Vi. a] Gii [{li-iy J- ^ . a] -E^ 27rif ^ 



-f^J Ea 3("«-9«)ra/3(n/3-g/3) 



(A.2) 



where we have discarded the A/ A self interaction and we defined the dual 
lattice charges Ua with a — i-\- (— |, |) (see fig. A.l) : 

^^.(Vj X Z) liy li_^^y 

and the quenched random dipoles 

q^ = i-V,.A = ^(Af - AU^ + A^- AU) 



27r 



27r 



Note that the neutrality of disorder charges Ea Qa = ^ follows directly from the 
definition of q (e.g. taking A to vanish at the boundary) and implies neutrality 
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Fig. A.l. Representation of the two networks. The plain and dashed hnes correspond 
respectively to the direct and dual networks. 

for the integer charges J2a na = (from the divergence of r(0)). The Coulomb 
potential T is defined on the lattice by 

A^r = r,- - 4Ti = v^Vo^r^.e, + Vj^v^^r^.e, = 27r 5^,0 (a.3) 

j=i±ex,i±ey 

which gives the expression 



i-j - -27r / ^2^^2 4 _ 2 cos(g^) - 2 cos(gy) ^ 



= r(0) + In |z - j\ + ln(2v^e^) + O 



(A.5) 



where the 7 is the Euler constant. Using the neutrality of the charges and 
the definition Ga-p = Fq — Fa-/3, and the definition of the disorder potential 
Va — —'^JY,pGap(l(3-i the lattice disordered CG action reads 

AcG = -pHcG = PJY1 ^aGapnp + /5 ^ n^Va (A.6) 



B Smooth cutoff procedure 



The Coulomb gas, or its equivalent Sine Gordon version, can also be renor- 
malized using smooth cutoff procedures. We will not give details here but refer 
the reader to [48] (see also [14]). Let us simply point out how, in that case, 
the full disorder Vr can indeed be decomposed, as in (6), into two bona fide 
disorders. 

In the case of a soft cutoff, the continuum approximation of the lattice Coulomb 
interaction reads (instead of (5)): 
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Gl"' = 2./(g,^(l--(..r,) (B.l) 

where 0(0) — 1 and we will choose ^{x) a positive monotonously decreasing 
function of x. One has the asymptotic large r behaviour G'f^ ~ Inr + C(0) + 
0(l/r) (see e.g. [14]). One also has: 



Gf = - 27rd/ / (^^^^(1 - cos(q • r)) + 0{df) (B.2) 

where a = ae'^''. Thus one can write the cutoff dependent decomposition Vr ~ 
V>M + and one gets upon increase of cutoff: 



and v^"'^ — v^"'^ + dvj. with: 

2j; f -ag0'(ag) ic 



dvrdvr' = Aiiardl / — ^' e^q-^ (B.4) 



Thus both T/^ and di'r are well defined physical gaussian disorders since their 
correlators have positive Fourier transform. In addition dVj. is short range 
correlated. For instance, taking (f){x) — e'^ one finds: 



^ ^rf/ e-i^-r'?K2a-) (B.5) 



C Renormalisation of the replicated Coulomb gas 



In this appendix, for completeness, we explicitly renormalize the vector Coulomb 
gas defined by (13). This amounts to extend to m-component vector charges 
the renormalisation of scalar CG [6]. The partition function reads: 

/ N j2 

^latt = E E n / ^>^[n.]e^-[---l (C.l) 

iV>2{ni...n;v}«=l ° 

where as usual the primed sum is over all distinct neutral charge configurations 
and the notation h.c. stand for all the hard core constraints |rj — r^l > oq for 
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all pairs implicit in the following. The action Aao is defined by 

AJni ...n^] = W 2X„, < In OiLZiA] n] (C.2) 

^ i^j \ do J 

A common way to renormalise usual Coulomb gas consists in coarse graining 
the partition function, leaving the expansion in number of fugacity (here Y) 
unchanged. This amounts to define scale dependent replica stiffness Kabil) 
and fugacities l^fn]. We will follow this scheme in this appendix. Note that 
another equivalent way would be to renormalise the correlation function di- 
rectly, following e.g [49]. As in the scalar case [2,6], renormalisation of the 
generalised vector Coulomb gas [50,51] proceeds in the same three steps : 
rescaling, fusion and annihilation (screening) of small dipoles. We now turn 
to the description of these three contributions : first we increase the hard-core 
cut-off Co — > Oo = OoC'^' with dl <^1. 



C. 1 Rescaling 

This increase of cut-off produces a naive rescaling : 

N j2 N j2 

TT / ^LZiy[nJe^-["i-"^l = ff / ^^lF[n,]e^^°["i-"^le'^'(2-^«''">') (C.3) 

where we used the neutrality J2i'>^i — ^ to express the correction coming from 
the action. We can absorb the extra factor in (C.3) to all order in Y by the 
change of fugacities corresponding to the equation 

diY[n] = (2 - Kabn''n')Y[n] (C.4) 



C.2 Annihilation (screening) and Fusion of charges 

Upon the increase of cut-off, two charges rip and have to be coarse grained 
if they are located in and with < |rp — r^l < agc'^^. Within the small 
charge density hypothesis, we consider only one such pair. For a dipole, these 
two charges have to be integrated out at scale do : this corresponds to the 
annihilation, while for a non neutral pair, the coarse grained charge is simply 
the sum of the two charges at scale Co (fusion). In both cases the partition 
function splits into Z = Z' + Zp^g where Zp^g ~ 0{dl) involves configurations 
with one pair of charges n^, distant of less than So while Z' doesn't. Zp^g 
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can be written as 



^v.- t n ^^y[^^] 



X E / 4j^4|^F[n^]y[n,]e-^-["-"-l (C.5) 

« J ar,<\r„—r„\<ar,e'^^ (l„ (l„ 



where, with the notation ccij = X^„ (,'^?2ir"^n^, the action reads 

(C.6) 

We must now distinguish between a neutral pair and a non neutral one. 
C.2.1 Fusion 

In this case the small pair of charges (np, riq) at scale Oo gives an effective 
charge Hp + at scale do-, located in R = (r^ + rg)/2. Thus we must integrate 
over the relative position of the two charges p = — when coarse graining 
the Coulomb Gas. To obtain the corresponding correction to Z' (of order one 
in dl), it is enough to expand (C.6) to order in p. This expansion reads 

{■n.i...Ta.f^} i=l,..N,i^p,q o "pjiq o \ "■o_k_"o'> q 



X 



Rl \ Oiin~\~OiiQ / I I \ OJj 



n ^ n . „ . 

Using ai^p + ai^g = ai^p^q, we can rewrite this correction to Z' as a single 
contribution to the fugacity Y[np + rig] of the non zero charge 11^ + 11^: 

diY[np + Uq] ^271 J2 ^W^K] (C.7) 



Note that in this expression rip and are distinguishable charges : this ex- 
plains the factor 2 between (C.7) and (14b). 



C.2.2 Annihilation 

In this case the small dipole of size p is integrated out at scale do- When coarse 
graining, we sum over both p and R, yielding a factor from the integration 
measure. Thus diverging contributions correspond to the expansion of (C.6) 
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to order 2 in p. Using ai^p — —ai^q, this expansion of Zp^q can be expressed as 



X 



{ni...njv} i=l,.-'V,i^p,<J o i<j¥=P,g 

I ^ / ^ I ^ + E «i,P«.-,P7P- V ln(r, - R)p. V ln(r,- - R) 



Performing the integral and reexponentiating the last term using the neutrality 
of the configuration {ni,...njv}, we get the correction to Z' coming from 
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c E n /^^N 



{ni...n;v}i=l,-Af,i7^P,g ° 



X n 



The constant C corrects the intensive free energy while the second term can 
be absorbed in a correction to the coupling constant 

diK~^^ = 27r2^nVF[n]F[-n] 



The three above contributions can be summarized into the set of RG equations 
given in the text (14) valid for all non zero vector charge, and for all m. The 
coefficient Ci and C2 depends on the IR regularisation. For our hard cut-off, 
we find within our procedure that ci = 2tt'^ and C2 = n. Note that the ratio 
ci/cg is independent of a uniform rescaling of the fugacities and is known, 
in the case of single component charges to be universal at a transition [49]. 
For a discussion of regularisation of replicated Coulomb Gas, see discussion in 
appendix F. 



D m — > limit of the replica RG equations 



In this Appendix we explicitly perfom the m — > limit of the whole set of 
RG equations (14a, 14b) with the restriction that the vector charges have only 
0, ±1 components. This hmit is taken using the parametrisation (16) of the 
fugacities Y[n] in terms of the function $(^+, z^). In the three different terms 
of (14a, 14b), corresponding to rescaling, annihilation and fusion contributions 
( see appendix C), we first modify sums to include fugacities for null charge, 
translate the expression in terms of $ and naturally take the m ^ limit. 
We will use the notation = / ^_ A ^{z+, z^). 
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D. 1 Rescaling 



The term corresponding to rescaling in (14b) is 



diY[n] = {2-naK'^W)Y[n] (D.l) 

which holds for any vector charge n. The second term can be expressed in terms 
of $ using (16) and K"-^ = K6"'^ — aK"^ with K = f3J . However the expression 
is much simpler if one uses, instead of (16), the equivalent parametrisation 
y [n] = e^P'^e^'^'"^{u, v) with p = n+ + U-, q = n+ - ri- and z± = e^("=^"). 
With l>, this yields 

Y^naK'^^'ub Y[n] = {(5Jp - a(5^J^q^) [ $(«, ^;)e'^f"e^«'' (D.2) 

a,b 



$(«, v)(j^- aJ^-^] e^f^e^^'' (D.3) 
\ OU Ov^ J 

^ _ / e'^P^e^'?" + (tJ^-^] H% v) (D.4) 

Ju,v \ OU OV^ J 

As this is true for any vector charge n satisfying = 0, ±1, i.e for any p, q or 
equivalently n+, n_, to satisfy (D.l) we can search for a function $ such that: 

In terms of ^{z+, Z-), this corresponds to the equation: 



di^z+, z_) = (2 + 0)^z+, z_) (D.5) 
= (2 + PJ{2 + z+d,^ + z.d,_) + <T(5''j\z^d,^ - z^d,_f) $(z+, z^) 

where we have used that $(«, v) = 2zj^z^^{zj^^ z_). Note that for this process, 
the integral of $, A/" = / $ satisfies simply diM = 2J\f. 



D.2 Fusion 



In (14b), the fusion term of two charges n' + n" = n is restricted to n' 7^ 

and n" 7^ 0. Furthermore since n = corresponds to the annihilation which 
is treated separately in (14a) (see below) and must not be counted twice, 
the equation (14b) can be used only for n 7^ 0. It is convenient to extend 
the equation (14b) to include Y[n = 0] = / $ = jV for which the fusion 
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contribution should be absent. Thus the extended equation corresponding to 
fusion reads: 



diY[n] = C2 E Y[n']Y[n"] - 5„,o E y[n']Y[-n' 

n'+n"=n n^O 
\ n',n'VO / 



(D.6) 



= cj E Y[n']Y[n'']-2Y[O]Y[n]-5^JT.y[^'n-^']-m0?]] 

\n'+n"=n \ n' / / 

where in the second equahty we allow for n' = or n" = 0. Turning to the 
representation in terms of n+, n_ we have: 



E Y[n']Y[n"] 

n'+n"=n 



Using permutation symmetry we find that (D.6) can be written as: 



diY[n] = 



z'_ + z'L 



l + z'_z'i + z'+z'L 



-2c2M{zl+zl- 



\1 + zLz'l + z'+z'l J 



The following choice for di^{z+, z^) allows to satisfy the above equation for 
all n (up to now m is still arbitrary): 



z'+zl 



z' +z" 



m{z+,z-)^ 

c. ((1 + z'.z'i + z'^z'^rs - ' ' ' ^ 

-2c2J\f^z+, Z-) - C2S{z+)S{z-) (((1 + z'_z'l + z'+z'^r - l)^^ - -A^') 



(D.7) 



l + z'_z'l + z'+z'lj/^^ 



We can now take the limit m ^ explicitly on this equation, which yields: 



di<!>{z+,z-)=C2(5 [z+- 



1 + z'_z'^ + z'+zl ^ 



5 U- - 



i + z'_z'; + z'+z'Lj/^^ 



-2c2J^^z+, z_) + C25(z+)5(z^)X^ 



(D.8) 
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D.3 Annihilation 



The screening equation (14a) for the rephca couphng constant reads 
which corresponds to the RG equations for a and J 



(D.9) 



di{TJ-^ + (7) = - ^ ^ n"n"F[n]r [-n] 



diG = 



Cl 



a=l n^O 

m 



m{m-l) „^i^o 



(D.IO) 
(D.ll) 



The sums over the charge n in (D.IO, D.ll) can be easily expressed in terms 
of $ using the variable n+,n_. The sum in (D.ll) thus reads 



J2 ^nV^"F[n]F[-n] 

a^b=l n^O 

E C;;,-- [(n^-n^r-(n^ + n_)]{(z'^z'Lr-(z'_z';r-) 

0<n-|_+n_ <m 

- z'_d,,y - + z'_d,>j 

X E t^;;+'"-(4^::)"M^-4)"-)) 

0<n++n_<m / ^(^2:')^{z") 

= m{m - 1) {{z'^z'L - z'_z'[ni + z'^z'L + ^'-4)""-'),^,^^^,,^ 

where C^+'"- = m\ / n^\n_\{m — n+ — n_) \ . Taking the m — > limit of the last 
equation yields directly with (D.ll) 



dia = Cl 



z\z" - z' zl 



1 + zLz'L + zLz'L 



(D.12) 



/ *(z')*{^") 

By the same method we get for the scaling equation of 



di{TJ 



Cl 



' z\z"_^ z'_zl^ ^z\z"_z'_zX 
(1 + z'^z'L + z'_z%Y 



(D.13) 



*(z')*(z") 
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D.4 Final set of RG equations 



Putting all these contributions together, we obtain the scaling equations for 
the coupling constant J, the correlated disorder strength a and the distribution 
of local disorder $(2;+, z^) : 



z'z" -z' z'[ ^ 



2_) = (2 + (D.14c) 
Li z'^z'L \l z'_ + z'L \\ 



where the diffusion operator has been defined in (D.6). 

E Connection between the direct and the replica method 



In this appendix we explore the connections between the expansion in number 
of sites of the free energy, and the expansion of the replicated partition function 
in power of composite charge fugacities Y[rL\. As we will see, both expansions 
coincide exactly, and we can thus consider the expansion of the replicated 
partition function as a generating functional of the site-expansions. 



E.l The 2 point free energy f^"^^ 



Fist we consider the first term /^^^ given by (30) of the expansion in number 

of independent sites (28). This term corresponds exactly to the approximation 
of independent dipoles considered in [23] : dipoles do not interact with each 
other, and the free energy is thus the sum over all positions of the pairs of the 
free energy of a pair : F — J^ry ffli ■, ^-^d the free energy of a dipole is simply 



= ln(l + Wry) ■ Wry = ( ^— ^ {Wry + W, 



r,r' I "^r',ry 

ao 



Wr>,r = e'^^^'-^r' 
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Using the decomposition of the disorder (section 2.2) Vr — + Vr, and 
averaging over the correlated disorder using rephca, we obtain 



rV> 7TT — ^ T-V> 



ln(l + Wry) ~m^o (1 + Wr,r')^ - 1 (E.l) 

/ Ir - r'l \ -2(P+9)/3-^+2{p-9)V/32j2 
= {z^Z^riz^Z^y l^—^] (E.2) 



l<p+q<m. 



where we used the non local correlation of the disorder (V^ — V^)"^ — AaJ'^ In (^^^-^ 
and the definition z^. = ye^^'"'' . With the bare definition of the replica charge 
fugacitics, this finally gives after average over Vr, the expected first term of 
the expansion of Z"* : 



ln(l + WryY' ^m^o E >^[n]r[-n] ( ) (E.3) 

[n^O] V °o / 



E.2 The 3 points free energy f^^^ 

The methods is the same as in previous Section, but the calculations are 
slightly more tedious. We consider the second term of the free energy expansion 

/ 1 + Wry + Wry + Wr'y \ , . 

V(l + Wr,r'){l + Wr,r")(l + W^r',r") J ^ ' ^ 

With the help of the previous Section, it is enough to consider only 

In (1 + Wry + Wry + Wr'y) 
Using the decomposition of figure (E.l), we obtain the sum 

In (1 + Wry + Wry + Wr'y) 

~^^0 (1 + Wry + Wry + Wr'yr - 1 



{p,q,s) pi-\-p2=P qi'\~q2=q si+S2=s 

0<p-\-q-\-s<m 



r — r \ / r — r \ / r — r 



Qjn I \ do I \ Ojfi 



By averaging the second term over , and using Zj^ — ye'^^'"", we get 



75 




Fig. E.l. possible configurations of neutral triplet of vector charges and parametri- 
sation of this triplet by the numbers of pairs of components (+1, —1): (p,q,s) 



v> 



2<7/32j2(pi-p2)2 



X 




z.z 



z z. 



\S2 



2a/32j2(gi-Q2)2 



|r'-r"| 



2ct/32j2(si_s2)2 



Inserting this result in the above expression averaged over V^^, the sum over 
the partitions of the interval [0, m] in < Pi +P2 + 9i + Q'2 + si + S2 < m can be 
exactly rewritten as a sum over all neutral triplets n, n', n" of m component 
vector charges with components 0, ±1 (see figure E.l). Thus one recovers the 
expression 



J2 Y[n]{r)Y[n']{r')Y[n"]{r")e-^",n',n" _ i 

n+n'+n"=0 



We now note that the last three logarithmic terms in (E.4) just give additional 
restrictions on this sum (as they correspond respectively to n = 0, n' = 0, n" = 
0). Thus, averaging over the local fugacities, we end up with the expected 
second term (see 11) of the expansion of in power of the vector fugacities 
y[n]: 



J2 Y[n]Y[n']Y[n"]e-'^"- 



n,n' ,n" 



(E.5) 



n,n' ,n' 



These simple combinatorics can be done on higher order terms : it gives the 
equivalence term by term between the m ^ limit of the expansion in ^[n] 
and the expansion in the number of independent sites of the moments of the 
free energy. 
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F Regular izat ion of replicated Coulomb gases 

In this Appendix we discuss the consequences on the RG equations of the 
choice of cutoff made in this paper for the rephcated Coulomb gas. We illus- 
trate for simplicity only the case of zero disorder. Although we will be mainly 
concerned with Coulomb gas, most of this discussion can be applied to other 
general replica field theory. 

Let us first recall the results for a single component CG (m = 1). We restrict 
to the most relevant charges ±1 for simplicity. It is defined by the action 

-(5H^K Y: n(r)G'(|r-r'|)n(r')+ln(y)^n^(r) (F.l) 

\r-T'\>ao 

The corresponding RG equation was derived by Kosterhtz: 

diy{l) = {2-K)y + 0{y') (F.2) 

We now consider m copies of this model in the absence of disorder, as illus- 
trated in Fig. (F.l). They are a priori physically completely uncoupled. The 
most natural cutoff procedure in that case would be independent cutoffs (e.g. 
hard core for each) (left figure). Another procedure, which becomes much more 
convenient in the presence of disorder, is to reformulate the m copies as a sin- 
gle Coulomb gas of vector charges with m components. However in that case 
the cutoff is by definition columnar (right figure) (e.g hard core vector charge 
are a hard columnar disk) and in a sense the copies are coupled, via the cutoff. 
We now check that in the pure case the ensuing vector CG RG equations are 
still perfectly compatible with (F.2) as they should. They read: 

9iF[n] = (2 — Kn ■ n) Y[n] + C2 ^ F[n']F[n — n'] -|- higher order terms 

Let us first illustrate the case of two copies m — 2. We can choose ^[1,0] = 
r[-l,0] = r[0,l] = Y[0,-1] = and = = = 

F[— 1, — 1] = Y2. Then considering all the possible fusions within this set the 
RG equations read: 

diY^ = {2-K)Y^ + Ac2Y^Y2 
diY2 = (2 - 2K)Y2 + 2C2YI 

A solution of these equations, to the order o(y^) at which we are working, is: 

yi[l,0]=Aiy, Yi\\\\=\2y] 



77 



with A2 = C2XI, where yi satisfies the single copy equation (F.2). This can be 
generahzed to higher charges so that in general one can find solutions of the 
type: 

= ^[ i,i, . ..,i ,o....o] = ^^y'l 

with coefficients Xp which can be determined for a given regularisation proce- 
dure. They depend on both the initial (m = 1) and the replica regularisation 
(see fig F.l). 




R R 



Fig. F.l. Schematic representation of the definition of the replicated vector charges 
and their cutoff 

This example illustrates how a generic cutoff procedure will produce non trivial 
coefficients Aj. Thus if one reinterprets formally Y[n] =< >$(y) in 

terms of a "disorder" $(?/) as done in this paper one must keep in mind 
that even the pure system corresponds to a non trivial "bare disorder" in the 
fugacities, which solely originates from the (convenient) choice of a columnar 
cutofiP^ It is thus clear that the definition of the local fugacity distribution 
is strongly cutoff dependent. 



G Explicit derivation of the expansion in the number of points 

In this appendix we derive formula (29) for the expansion of the free energy 
in the number of points. The same method can be applied to other physical 
observables. 

To organize this expansion we start by introducing fictitious site dependent 
fugacities for the charges : Crr These are introduced only as a trick in this 
Appendix and should not be confused with the real disordered fugacities of 
(6). Indeed, to recover the original model, we will set them back at the end 
either to Cr, = 1 (for the lattice model (4)) or = y for the continuum model, 
where y is the corresponding fugacity for the charge in the pure case. These 



with = Z- = y 

in the interpretation as a branching process, it does correspond to the disorder 
in the tree structure. See the end of Section 5.5 for a particular choice which does 
not introduce additional disorder 



78 



fugacities are introduced by writing a more general form of the the partition 
function (4) as 

Z[VX] = Y. t E (nCf)e^^^^»-."'^--^"^^^^'"'^' (G.l) 

P {ni,...np} ri^---5^rp \ i / 

Here and below, as in (7), all formulaes can be extended to the continuum 
model by replacing discrete sums over distincts sites ri ^ v^ - ■ ■ 7^ by 
integrals with, e.g hard core conditions |ri — > Co- Note that in the above 
expression (G.l) we do not make use of the decomposition (6) and V denotes 
the original disorder. 

Let us consider for simplicity a system of distinct points ri, ...rjv. The free 
energy functional C] — —2^ In C\ is a function of the N variables Cri,--) 
Criv • Let us write the conventional Taylor expansion of the free energy around 
C = 0: 



P1,P2,..PJV=0 ■■■PN- ^^^^ ■ ■ ■ ^^^M 

We now separate in this sum the terms which have only one non zero pi, then 
only two non zero pi, etc.. Thus we can rewrite: 



+00 

= E E CJ^C] 

where the sum is over all distincts sets (r^j, .., Tj^} of A; distincts points (among 
the TV points) and the sum over each qt goes from 1 to +00. 

We have thus obtained an expansion of F[V, Q as a sum of terms of the form 
/rf rfe which depends exactly and only on the variables (and thus also only 
on the variables V) evaluated at the k distinct points r . . . r j^, . For a neutral 
Coulomb gas it starts with k — 2 and reads: 



mQ= E f^U[V,Q+ E fSU,r.,\V,C\ + ... (G.4) 
with the definition 



ft-r..\V,(]= E a.^dC^^ lc^oxC^...Q (G.5) 
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Note that in this last expression we can drop oTit the dependence of F[V, C] 
on the fugacities (r and the potential K- at points different from r , . . . r . 

A more explicit expression can be obtained by summing over the qi in (G.5): 

/S.rJ^,C] = E(-l)'-' E ^r,,...r,[l-,C] (G.6) 

1=0 h,...iie[i,...k] 

where -^rij,...rijl^, C] is the free energy associated with the system of sites 
rj^,...rj, (instead of the full lattice). Equivalently, it does depend only on 
the fugacities Cr (or potential Vr) at points rj^,...rj;. After setting = 1 
(lattice model) or (^^ = y (continuum model) this gives the definition of the 
that we use throughout this paper, namely the equation (29). This 
last expression (G.6) allows to exphcitly compute the coefficient /^'^^[F] of the 
expansion (G.4) for arbitrary order k. In the Section (4.2) we use it to give 
explicit expressions for the first few terms of the expansion of the free energy. 



H Fusion corrections to the free energy expansion 

As shown in the previous Section the free energy expansion involves the par- 
tition function 2^} ^.^ of a system of finite number of sites. We illustrate here 
how the rule of fusion of environments described in Section 4.3 works more 
generally. Let us consider a system of p > 4 sites, with charges restricted to 
0, ±1 and of energy and fugacities given by (7). Upon increase of the cutoff one 
must take into account fusion of the sites ri and r2 into the site f = |(ri-|-r2), 
one finds that: 

- Z^.r, (H.l) 

with: 

+ E (^a„..r,+^it,..rj (H.2) 

q=2,q eveni2,...iqe[3,..p] 

where Wu = + ^t^^ri 

E E e-'''"^'^' (H-3) 

ni=±l nri=±l 
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where H has been defined in (7) and the charge ni is located in r. Note that 
w\^},^ = Wr^Y-i the dipole weight defined in Section 4.3.2. In (H.2) the first 
term corresponds to the total weight of configurations with no charge in ri 
and r2. The second term correspond to configurations with a dipole in (ri, r2). 
Its original expression is complicated but simplifies when ri and V2 are fused 
into f ( the interaction energy of any given other charge Ua in with rir^ 
becomes opposite to the one with rir^ (up to higher order terms) and Wr^,r2 
simplifies into W12 to lowest order in dl as explained in Section 4.3). The last 
two terms correspond to all configurations with one charge either in ri or in 
r2- 

The important property with respect to the free energy expansion is that one 
can factor the term 1 + W12 and rewrite: 

ln^£;^l, = ln(l + iyi2)+ln(Zg-)+ E ^S„..rJ 

q=2,q eveni2,...iqe[3,..p] 

where W^'^^ as the same definition as W^'^^ except that the fugacity at f has 
been modified according to the fusion rule for fugacities (42). The total factor 
inside the last logarithm is exactly the partition function Z^^^\^ of a system 
of p — 1 sites with the new fugacity given by (42) on the site f . 



I Higher charges and extensions 

In this Appendix we briefiy indicate how higher charges can be included in the 
same (fusion-diffusion) formalism and why including them does not affect any 
of our results. Let us consider for instance the charges ±2 and define n+_|_ and 

n respectively as the number of component charges +2 and —2 in the vector 

charge n. The parametrization (16) can be readily extended to encode also for 
the charges ±2 simply by writing the vector fugacity Y\p\ as an average over 
a function $(^+, 2;++, z ) as 

n_, = / $(z+, z++, Z--)z++ zl' z++^ zls 

J Z+,Z-,Z++,Z 

where the random variables 2;++ and z represent the random local fugac- 
ities for the charges ±2. Similar manipulations as in Appendix (D) lead to 

a RG equation for a normalized P(2;+, 2;++, 2; ) which we will not write 

explicitly. It does contain a diffusion operator as well as a fusion term. For 
illustration we simply give the diffusion operator, expressed using more con- 
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venient variables z± — e^^'^^'"^ and z++ = z^e^^"^' ^'"'^ and z = 2;^ e'^*^"' It 

reads, for the corresponding probability distribution P{u,u' ,v,v'): 



2-J{du + 2du,) + aJ^d^, (I.l) 

The detailed study of the corresponding RG equation, together with the fusion 
term will not be reported here. Instead let us indicate how the irrelevance of 
the higher charges can be justified. Since one has: 



2 - n.K.n =2- pj{n+ + n_ + 4(n++ + n__)) (1.2) 
+ap^J^{n+ -n- + 2(n++ - n__))^ 

it is clear that the reduced probability 



Q{z++,z—)^ P{z+,z-,z++,z—) (1.3) 

JZ+,Z- 

satisfies the same diffusion equation as P(2;_|_,2;_) in (21) but with the change 
J ^ AJ and a — > (j/4. It becomes relevant by power counting only at cr = 
4(7c — 1/2. Thus it is less relevant than P{z+, zJ) in the region of the phase 
diagram of interest and we can rightly neglect the new fugacities z^^ and 

z One can also check that the fusion terms imply that -P(-2++ ~ 1) is of 

order P(-2+ ~ 1)^ when this parameter is small (since at lowest order the 
RG equation contains a term proportional to ^(2;++ — z'j^z'\^)). The rare events 
which involve the charges ±2 are thus subdominant. 



J Evaluation of screening integrals at criticality 



To evaluate the integral (105a) in the screening equation for J we consider 
separately three intervals u — Xi < b, Xi + b < u < —Xi — b and u > 
—Xi — b, where 6 is a number such that pi{x) can be replaced by its asymptotic 
expression for x > b with good accuracy. Only in the middle interval can we 
use the universal tail expressions for both factors pi. Thus we have, using 
symetries: 



dupi{u - Xi)pi{-u - Xi) = 2 dvpi{v)pi{-v - 2Xi) 

J— 00 

+ r^^^'piiu - Xi)pi{-u - Xi) (j.i) 

JXi+h 
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In the first integral we can replace the second function by its asymptotic form 
and get: 



^g2X,/VS r dvpi{v){-2Xi - w)e^e-^T^ ^ CX,e'^' (J.2) 

since, remembering that the integral dvpi{v) — 1, the above integral is at 
most a finite number. The second integral in (J.l) is estimated as: 

and gives the leading contribution. For X;// — > it behaves as ~ Xfe^^'/"^ 
and gives the estimate in the text. It is encouraging to note that the leading 
behaviour precisely originates from the interval where pi{u) can be replaced 
everywhere by its universal asymptotic form. 

The integral (105b) in the screening equation for a can also be estimated by 
considering for each variable u, u' two intervals u < b and u > b. The end 
result is that we find to leading order: 

diJ ~ -^dia ~ A(_X;)3e2^'/v^ dv{l - ^;2)e-S(i+-') (J.3) 
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